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Abstract 

After a short survey of the physics of solar neutrinos, giving an overview of 
hydrogen burning reactions, predictions of standard solar models and results 
of solar neutrino experiments, we discuss the solar-model-independent indica- 
tions in favour of non-standard neutrino properties. The experimental results 
look to be in contradiction with each other, even disregarding some experi- 
ment: unless electron neutrinos disappear in their trip from the sun to the 
earth, the fluxes of intermediate energy neutrinos (those from 7 Be electron 
capture and from the CNO cycle) result to be unphysically negative, or any- 
way extremely reduced with respect to standard solar model predictions. Next 
we review extensively non-standard solar models built as attempts to solve 
the solar neutrino puzzle. The dependence of the central solar temperature 
on chemical composition, opacity, age and on the values of the astrophysical 
S-factors for hydrogen-burning reactions is carefully investigated. Also, pos- 
sible modifications of the branching among the various pp-chains in view of 
nuclear physics uncertainties are examined. Assuming standard neutrinos, 
all solar models examined fail in reconciling theory with experiments, even 
when the physical and chemical inputs are radically changed with respect to 
present knowledge and even if some of the experimental results are discarded. 



1 



Typeset using REVTeX 



Contents 



Introduction 



Uverview| 



A Hydrogen burning reactions! 



B Stellar structures and standard solar models 





c 


Results of standard solar models 




L) 


Solar neutrino spectrum and the predicted neutrino signals 


E 


A look at experimental results 




F 


A remark on errors 








11 


JNeutrino Muxes almost independently of solar models 




A 


The solar luminosity constraint 


B 


The minimal Gallium signal for standard neutrinos 




C 


Chlorine and K AMiOK AN DL 



D The space for intermediate energy neutrinos 
E Four equations and four unknowns 



F Where are Be and CNO neutrinos? 



G Probabilities, confidence levels and all that 



1 Neutrino mass and (3 decay 



|2 Statistics and the solar neutrino problem 



H Experimental results and standard solar model; 



I The beryllium-boron relationship 



J What if the sun were burning less no"wi 



K Universal neutrino oscillations] 



L Concluding remarks 



III Non-standard solar models: why and how? 



A Introduction 



B Cooler solar models: why and how'; 



C Homology relationships 



D Scaling laws for the physical quantities 



E Scaling laws for neutrino fluxes] 



Dependence on the central temperature 



IV Low central temperature solar models 



A Introduction] 



B The p+p — > d + e + + v reaction rate 



C Radiative opacity 



D" The "metallicity" of the sun 



E Rejuvenated suns| 



F Low temperature models and neutrino fluxes: summary and explanation 



C More about homology: additional consequences and a possible test 



H On the accuracy of the central solar temperature 



4 

7 

8 
9 

12 
14 
16 
18 

20 

21 

22 
23 
24 
26 
28 
28 
29 
31 
33 
33 
34 
35 
36 

37 

37 
37 
39 
41 
42 
43 

44 

44 
45 
46 
47 
48 
49 
51 
52 



2 



V Nuclear reactions in the sun 



53 



A 


The status of d He +^He - 


> 4 He + 2p reaction 


B 


'the status of °He + 4 He - 


-> 'Be + 7 reaction 



c 



55 

Neutrino fluxes and helium reactions! 55 

58 

59 

60 

60 



D 


The p+''Bc - 


-> 8 B +7 reaction 


E 


The p+ i4 N - 


-^O +7 reaction 



Plasma screening of nuclear charges 



1 The results of different models 



2 A model independent analysis! 62 



VI JN on-standard solar models and experimental results 



64 

Introduction! 64 

65 

66 

66 



B Fitting all experimental results 



C What if one experiment were wrong 



D What went wrong with solar models' 



Vll The never ending story? 
[APPENDIXES 



~K Our standard solar model 



B The stability of pep and CNO neutrino fluxes 



References 



68 
69 
69 
71 
72 



3 



INTRODUCTION 



In a prophetical paper of 1946, Bruno Pontecorvo [T|] reviewed the arguments for the 
existence of neutrinos (at that time based mainly on conservation laws) and stated: 

"Direct proof of the existence of the neutrino . . . must be based on experiments, the in- 
terpretation of which does not require the law of conservation of energy, i.e., on experiments 
in which some characteristic process produced by free neutrinos . . .is observed." 



Then he pointed out the relevance of inverse beta decay for neutrino detection and 
outlined the physics potential of neutrinos from nuclear reactors. He also mentioned solar 
neutrinos: "The neutrino flux from the sun is of the order of 10 w neutrinos/cm 2 /sec. The 
neutrinos emitted by the sun, however, are not very energetic . . . ". In 1946 Pontecorvo was 
not optimistic about solar neutrino detection. Since then, solar neutrino physics has made 
enormous progress, with a significant acceleration during the last decade. 

From the seventies to the first half of the eighties the only result was that of the Chlorine 
experiment 0, essentially implying an observed 8 B neutrino flux lower than that predicted 
by Standard Solar Models (SSMs). Many physicists were not convinced of the existence of 
a solar neutrino problem, mainly because of the uncertainties in estimating the flux of the 
rare, high energy 8 B neutrinos. In 1989 || results from the KAMIOKANDE experiment 
were presented, confirming the low 8 B neutrino flux. However, as observed by Bahcall and 
Bethe , the comparison between KAMIOKANDE and Chlorine data implied an additional 
puzzle. As the Chlorine experiment is also sensitive to 7 Be neutrinos, if one subtracts 
from the CI signal the 8 B contribution as implied by KAMIOKANDE data (for standard 
neutrinos), no room is left for 7 Be neutrinos. Thus the problem started involving also these 
less energetic neutrinos, for which the theoretical predictions are much more robust. At the 
beginning of the nineties, the results of Gallium experiments (GALLEX || and SAGE ||) 
became available; for these the signal is expected to depend primarily on pp neutrinos, but 
to a lesser extent also on 7 Be neutrinos and on 8 B neutrinos. In both Gallium experiments 
the signals were again significantly smaller than the SSM predictions. The exposure of the 
Gallium detectors to a 51 Cr neutrino source |7]|| was important for confirming quantitatively 
the sensitivity of the apparatus for neutrino detection. 

By combining the results of the various experiments, it became possible to obtain relevant 
information on neutrinos independently of solar models. As we shall see, the solar neutrino 
experiments tell us that the following assumptions are most likely contradictory: 

• Neutrino properties are correctly described by the minimal standard model of elec- 
troweak interactions (i.e. "standard" neutrinos). 

• Solar energy arriving onto earth is produced in the sun by nuclear reactions of the 
form: 

Ap + 2e~ ^ 4 He +2u e . 

Since there is little doubt that the sun, as many other stars, is powered by the conver- 
sion of hydrogen into helium, then experiments point towards some non-standard neutrino 
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properties and — as was the case in 1946 — energy conservation becomes again a key tool, 
albeit an indirect one, for the study of neutrinos. 

This paper attempts to answer the following questions: 

i) which information on neutrino properties can be obtained directly from solar neutrino 
experiments, independently — as much as possible — of solar models? 

ii) how much room — if any — is still open for an astrophysical or nuclear physics 
solution of the solar neutrino puzzle? 

As the SSM has to be abandoned in face of the clear disagreement with all experiments, 
we shall present a systematic analysis of non-standard solar models (i.e. models where 
some input parameter is varied beyond its estimated uncertainty), in order to explore the 
possibility of reconciling theory and experiments. 

As we are writing this paper, a generation of solar neutrino experiments has 
been essentially completed. In preparation for the exciting future experiments (SU- 
PERKAMIOKANDE [|, BOREXINO (TTJ, SNO Q, HELLAZ JT2), ...), it is time to 
summarize what we have learnt so far. 

In the minimal standard model of electroweak interactions, neutrinos are massless and 
consequently there is no mixing among weak flavour states. Also their magnetic moments 
vanish and they are stable. The experimental hint of non-standard neutrino properties 
(i.e. masses, mixing, magnetic moments, possible decay schemes, . . . ) thus seems to indicate 
some physics beyond the standard model. The present situation looks to us very similar to 
the one described by Pontecorvo in 1946: just change the word existence to non-standard 
properties, in the sense that solar neutrino experiments together with energy conservation 
strongly point towards non-standard neutrino properties. Much as in 1946, the future of 
neutrino physics now demands experiments yielding decisive direct evidence of non-standard 
neutrino properties. 



The plan of this paper 



In the first section we present a short and simple introduction to the field for non experts, 
reviewing: i) hydrogen burning in the sun and neutrino production, ii) the predictions of 
Standard Solar Models (SSM) and iii) the results of solar neutrino experiments. 

In the next section we discuss the information on neutrino fluxes which can be derived 
directly from experiments, almost independently of solar models, assuming standard neutri- 
nos. This is the main part of our review, as it is meant to demonstrate that, for standard 
neutrinos : 

• the available experimental results appear to be mutually inconsistent, even dismissing 
one out of the four experiments. 

• Even neglecting these inconsistencies, the flux of intermediate energy neutrinos 
(Be+CNO) as derived from experiments is significantly smaller than the prediction 
ofSSM's. 

• The different reductions of the 7 Be and 8 B neutrino fluxes with respect to the SSM 
predictions are essentially in contradiction with the fact that both 7 Be and 8 B neutrinos 
originate from the same parent 7 Be nuclei. 
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In Sec. f f| we describe a strategy for building non-standard solar models leading to 
reduced 7 Be and 8 B neutrino fluxes. There are essentially two ways: i) producing models 
with smaller central solar temperatures; ii) playing with the nuclear cross sections which 
determine the branches of the fusion chain. In the same section we introduce some algorithms 
for discussing the non-standard models, to be presented in the subsequent sections: we 
essentially seek parameterizations of the results (neutrino fluxes and physical quantities 
characterizing the solar interior, . . . ) in terms of the input quantity which is being varied. 



In section |TV] we examine non-standard solar models with a central temperature T c 
different from the SSM prediction. Clearly the central temperature is not a free parameter, 
its value being fixed by the equations of stellar structure and by the physical inputs one is 
relying on. We present a systematic analysis for the variations of several physical inputs 
(S pp , opacity tables, chemical composition, age, . . . ) affecting T c . An interesting feature 
common to the large majority of the models that we examined, is that the radial temperature 
profiles appear largely unchanged, aside from a scaling factor which can be evaluated, e.g., 
as the ratio of central temperatures. Thus, the profile is the same for a given value of T c 
independently of which parameter is varied in order to obtain that central temperature. 
This allows one to compare experimental data with the theoretical predictions as a function 
of just one parameter, namely T c . 

Section is devoted to the role of the nuclear cross sections which can influence the 
branching ratios of the p-p chain. We discuss in detail the role of the He+He reactions and 
we consider possible effects of to screening of nuclear charges by the stellar plasma. We also 
discuss the role of the p+ 7 Be reaction. 

In section |VT| we compare the non-standard solar models just built with experimental 
data. We also consider hybrid models, where both T c and nuclear cross sections are varied 
simultaneously. All our attempts to reconcile theory and experiments having failed, we 
conclude that 

• for standard neutrinos, the present experimental data look in disagreement not only 
with the Standard Solar Model, but with any solar model that we are able to build. 

In the Appendix A, we specify our Standard Solar Model in some detail. We concentrate 
here on our own calculations as we obviously have for these a more detailed knowledge 
of the physical and chemical inputs and of outputs of the code. In fact, several groups 
have produced in the last few years updated and accurate standard solar model calculations 
|T3^p3|. In this paper we shall not give a systematic presentation of these. We also do not 
describe solar neutrino experiments nor helioseismology. All these matters are excellently 
covered in other review works, e.g., 

- Bahcall's book [24j] and papers in Review of Modern Physics fll3-15fl, all extremely useful, 



centered around standard solar models. 



- The report by Koshiba ||25|| , providing a fascinating journey through the experimental 
techniques of neutrino astrophysics. 

- The review by Turck-Chieze et al. [p6| , again centered around standard solar models, 



including a clear discussion of helioseismology in relationship to the solar neutrino problem. 
- A recent reprint collection p7 |, presenting most of the significant papers on solar neutrinos, 



and an excellent bibliography. 
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I. OVERVIEW 



The discovery of the enormous energy stored in nuclei early led astrophysicists to spec- 
ulate that reactions among nuclear species were the source of the energy in stars: 

" Certain physical investigations in the past year make it probable to my mind that some 
portion of sub-atomic energy is actually being set free in a star. . . . Aston has . . . shown 
. . . that the mass of the helium atom is less than the sum of the masses of the four hydrogen 
atoms which enter into it. ... Now mass cannot be annihilated, and the deficit can only 
represent the energy set free in the transmutation. . . . If five per cent of a star's mass con- 
sists initially of hydrogen atoms, which are gradually being combined to form more complex 
elements, the total heat liberated will more than suffice for our demands, and we need look 
no further for the source of a star's energy . . . (Eddington 1920) p8|P|. 

After the discovery of the tunnel effect in 1928 p9|j30||, Atkinson and Houtermans pre- 



sented the first qualitative theoretical approach to the problem |yj . From the nuclear data 
available to them and the already known fact that hydrogen is the most abundant element 
in the sun and the universe, they concluded that the energy-producing mechanism involved 
primarily hydrogen, as suggested by Eddington. As a matter of fact, the fusion of hydrogen 
nuclei to form helium is still the only known process that can supply the required power for 
the long solar life. 

Hydrogen burning can be represented symbolically by the "fusion" reaction: 

Ap + 2e- -> 4 He + 2u e . (1) 

The total energy released in Eq. ([!]) is Q= 26.73 MeV and only a small part of it (about 
0.6 MeV) is carried away by the two neutrinos. 

From the solar radiative flux at the earth, called the solar constant K & = 8.533(1±0.004)- 
10 11 MeV cm -2 s _1 |15|], one immediately derives the total flux &tot °f electron neutrinos 
arriving on the earth, if they are not lost en route (by neglecting the small fraction of energy 
carried by neutrinos): 

$ tot ^ 2K G /Q « 7 • 10 10 cm- 2 s" 1 . (2) 

The definite proof of nuclear energy production in the solar core lies in the detection of 
solar neutrinos. Since the cross sections for neutrino detection depend strongly on neutrino 
energy, the energy spectrum of solar neutrinos has to be known and that requires a detailed 
knowledge of the reactions summarized by Eq. ([!]). 



In the same paper, a few line ahead, one reads: "If indeed the sub-atomic energy in the stars is 
being freely used to maintain their great furnaces, it seems to bring a little nearer to fulfilment our 
dream of controlling this latent power for the well-being of the human race — or for its suicide." 
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A. Hydrogen burning reactions 



In stellar interiors, nuclear interactions between charged particles (nuclei) occur at col- 
lision energies E well below the height of the Coulomb barrier and are only possible due to 
the tunnel effect, with a characteristic probability p9 |: 



exp 



(3) 



where Z^j are the atomic numbers of the colliding nuclei and /i^ is their reduced mass. The 
lightest elements thus react most easily and the conversion of H into He is the first chain of 
nuclear reactions which can halt the contraction of a new-born stellar structure and settle 
stars in the so-called Main Sequence phase. 

The mechanisms for hydrogen burning were first elucidated in the late 1930s indepen- 
dently by von Weizsacker and Bethe and Critchfield f3*2|-3%|; from their works it became 



clear that two different sets of reactions could convert sufficient hydrogen into helium to 
provide the energy needed for a star's luminosity, namely the proton-proton (p-p) chain 



and the CNO bi-cycle, which we are going to review briefly (see |35| for a more extended 
presentation). 

At the temperature and density characteristic of the solar interior, hydrogen burns with 
the largest probability (98%) through the pp chain. This develops through the nuclear 
reactions presented in Fig. [l|, where we show the different possible branches. 

The chain starts with a weak interaction process, either the pp or the pep reaction. The 
mechanisms for these two reactions are very similar and hence the ratio of pep to pp reaction 
probabilities is fixed almost independently of the details of the solar model; at the densities 
of the solar interior (p ^150 g/cm 3 ) the three body process has a smaller probability. The 
deuterons produced quickly burn into 3 He as a consequence of the much larger cross section 
of the electromagnetic process d+p— > 3 He +7 . 

After 3 He production, several branches occur. The hep reaction turns out to have a neg- 
ligible rate and will not be discussed further. Even if the 3 He abundance remains a factor 
10 4 lower than that of 4 He in the energy production region, the 3 He+ 3 He — ► 4 He +2p reac- 
tion is favoured with respect to 3 He+ 4 He— > 7 Be+7, because the probability of the reaction 
mediated by strong interaction is four orders of magnitude larger than that of the electro- 
magnetic process and because the tunnelling probability through the Coulomb barrier for 
the case of the lighter 3 He nuclei is a factor ten larger. 

With the 3 He+ 3 He reaction the chain reaches one of the possible terminations (pp-I 
chain). On the other hand, after the 3 He+ 4 He reaction, the chain branches again at the 
7 Be level, due to the competition between electron and proton capture. Note that, at 
the temperatures of the solar interior, the electron capture (although it is a weak process) 
dominates over the electromagnetic reaction since it has no Coulomb barrier. The chain 
passing through the 7 Be+e~ reaction is called pp-II chain, and the one involving the 7 Be+p 
reaction is called pp-III chain. 

To conclude this overview of the pp chain we remark that the basic process in any case is 
the aggregation of four protons to form a 4 He nucleus, the isotopes 2 H, 3 He, 7 Be, 7 Li and 8 B 
being intermediate products which are created and destroyed along the chain, their number 
densities staying constant (and small) when the chain reaches equilibrium. 



8 



Regarding the neutrino production we note that: 

i) in the pp-I chain only pp or pep neutrinos are produced, the latter with a small probability 
(0.2% in the standard solar model). 

ii) In the pp-II chain one pp (or pep) and one 7 Be neutrino are emitted. 

iii) In the pp-III chain one pp (or pep) and one 8 B neutrino are emitted. 

Figure |2] shows the main reactions in the CN and NO cycles which become efficient at 
rather high temperatures. The overall conversion of four protons to form a 4 He nucleus 
is achieved with the aid of C, N and O nuclei. The total energy release is clearly the 
same as in the pp chain. The ratio of the astrophysical S-factorsQ for the two reactions 
15 N (p, 4 He ) 12 C and 15 N (p,7) 16 being about a factor 1000, the role of NO cycle is generally 
marginal. In the sun essentially only the CN cycle (left side of Fig. ||]) is contributing to the 
energy production. 

For each fusion of four protons into 4 He via the CN cycle one neutrino from the decay 
of 13 N and another from that of 15 O are produced. Note that the efficiency of the process 
is determined by the reaction with the smallest cross section, which clearly is p+ 14 N — > 15 
+7, an electromagnetic process with the largest Coulomb repulsion. 



B. Stellar structures and standard solar models 

According to the theory of stellar structure and evolution, the condition for hydrostatic 
equilibrium jointly with the conservation of energy and the mechanism for energy transport 
determine the physical structure of a star and its evolution from four main inputs, namely: 

i) The initial chemical composition. 

ii) The equation of state for stellar matter. 

iii) The radiative opacity, k, as a function of density p, temperature T and chemical com- 
position. 

iv) The energy production per unit mass and time e, again as a function of p, T and chemical 
composition. 

For the equation of state (EOS) one has to evaluate the ionization degree and the pop- 
ulation of excited states for all nuclear species. In addition, one has to take into account 
several physical effects of the stellar plasma (like Coulomb interaction and/or electron de- 
generacy) which introduce deviations from the "perfect gas law" prediction. The study of 
EOS has improved along many years and accurate tabulations are available (see |I6| for a 
recent discussion). 

The radiative opacity k is directly connected with the photon mean free path, A = 1/ Kp. 
All throughout the internal radiative region k governs the temperature gradient through the 
well known relation [37|: 

dr 4acT 3 V ; 



2 We denote as Oij the cross section for the reaction between nuclei with atomic mass numbers 
i and j, is the relative velocity, E is the collision energy and the astrophysical S-factor is: 
Sij = <Jijexp[27rZiZje 2 /(hvij)]E. 
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where F is the electromagnetic energy flux. The evaluation of k as a function of the com- 
position of the gas and of its physical conditions requires a detailed knowledge of all the 
processes important for radiative flow (elastic and inelastic scattering, absorption and emis- 
sion, inverse bremsstrahlung, . . . ) and in turn a detailed evaluation of the atomic levels in 
the solar interior. As a consequence, the evaluation of k is a rather complex task, but the 
results have been continuously improving over the years. At present one relies mainly on 
the calculations of the OPAL group at the Lawrence Livermore National Laboratories (see, 
e.g. Refs. ||-|D|). 

Expressions for (the nuclear contribution) to e are derived essentially from the tables of 
nuclear reaction rates, originally compiled and continuously updated by Fowler, until his 
death (see Ref. jy] and references therein). 

According to Stix "The standard model of the sun could be defined as the model which 
is based on the most plausible assumptions." This means that physical and chemical inputs 
are chosen at the central values of experimental/observational/theoretical results. We prefer 
to call such a model the Reference Solar Model (RSM); throughout this review we use as such 
the "best model with diffusion" of Bahcall and Pinsonneault j0| denoted as BP95. More 
generally, let us define a Standard Solar Model (SSM) as one which reproduces, 
within uncertainties, the observed properties of the sun, by adopting a set of 
physical and chemical inputs chosen within the range of their uncertainties. 

The actual sun has a mass M = (1.98892 ± 0.00025) • 10 33 g @ and a radius R = 
(6.9598 ± 0.0007) • 10 10 cm ||44j| . Denoting by X, Y and Z the relative mass abundances of 



H, He and heavier elements respectively (X+Y+Z = 1), the present ratio of heavy elements 
(metallicity) to hydrogen in the photosphere is (Z/X) phot0 = 0.0245(1±0.061) J45],|T5|. At the 
age t© = (4.57±0.01) Gyr 0, the sun is producing a luminosity L = 3.844(1±0.004) • 10 33 
erg/sec [|T5|]. 

If a complete information about the (initial) photospheric composition was available and 
if the theory of stellar models was capable to predict stellar radii firmly, then there would 
be no free parameter. The theoretical model should account for the solar luminosity and 
radius at the solar age, without tuning any parameters. 

The photospheric helium abundance is however not strongly constrained by direct ob- 
servations, since no helium line exists at the photospheric depth, see Ref. |^6| for a critical 
discussion (helioseismology can however provide some indirect information, see later). Fur- 
thermore, the present observed photospheric composition is different from the initial one due 
to diffusion and gravitational settling. Moreover, the radii (but only the radii) of stars with 
convective envelopes (e.g. the sun) depend on the assumption made about the convective 
transport, as we shall discuss later on. As a consequence, one has the freedom of adjusting 
some parameters in the course of the calculation. 

In order to produce a standard solar model, one studies the evolution of an initially 
homogeneous solar mass up to the sun age. To obtain L , R Q and (Z/X) photo at t , one 
can tune three parameters: the initial helium abundance Y, the metal abundance Z and a 
third quantity related to the efficiency of convection (see below). 

The effects of these parameters can be understood simply. The luminosity of the sun 
(more generally of any star in the main sequence) depends in a rather sensitive way on the 
initial helium content Y; increasing it, the initial sun is brighter and a given luminosity is 
reached in a shorter time. Since the ratio Z/X is constrained by observational data, Y and 
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Z cannot be chosen independently: if Y increases, Z must decrease. 

To get the proper radius R , one adjusts the efficiency of the external convection, which 
dominates the energy transport in the outer layers of the sun (about the outer 30% of the 
solar radius). The precise description of the convection in the external part of the sun 
is an essentially unsolved problem, and the process is commonly described in terms of a 



phenomenological model, the so-called "mixing length theory", see e.g. Ref. [37|- Following 
this approach, we define the mixing length I as the distance over which a moving unit of 
gas can be identified before it mixes appreciably. This length is related to the pressure scale 
height H p —1/ (cflnP/dr) through 

I = a Hp, (5) 

where a is assumed to be independent of the radial coordinate and it is used as free param- 
eter. By varying a, one varies the mixing length and thus the convection efficiency, which 
determines the solar radius. Thus if a is increased, convection becomes more efficient, the 
temperature gradient smaller and the surface temperature higher. Since the solar luminosity 
is fixed, the radius has to decrease. 

In the last few years, one of the principal improvements of SSMs was the treatment 
of element diffusion g7^,|g^,0§I],|l3],^],|3|,|2| . Noerdlinger (1977) first included the 



effects of He diffusion on the evolution of the sun, while Cox, Guzic and Kidman (1989) 
were the first to take into account the settling of heavy elements. The stronger pull of 
gravity on helium and heavier elements causes them to diffuse slowly downward, relative to 
hydrogen. As a result, if diffusion is taken into account, helium and heavier elements in the 
solar photosphere are depleted with respect to their abundance in the original mixture. 

We remark that, as a strong simplification, the sun is taken as spherically symmetric, 
i.e. all physical quantities vary only radially. This relies on the assumption that the internal 
rotation is sufficiently slow [53] and the internal magnetic fields are sufficiently weak so that 



the corresponding forces are negligible. 

All in all, it looks that a solar model has three (essentially) free parameters, a, Y, and 
(Z/X)i n , and produces three numbers that can be directly measured: the present radius, 
luminosity and heavy element content of the photosphere. This may not look as a big 
accomplishment. At this stage, one's confidence in the Standard Solar Models actually rests 
on the success of stellar evolution theory to describe many, and more complex, evolutionary 
phases in good agreement with observational data. 

In recent years however, helioseismology has added important data on the solar structure, 
which provide severe constraints and tests of standard solar model calculations. Helioseis- 
mology can accurately determine the depth of the convective zone and the speed of sound 
0, at the transition radius Rb between convective and radiative transport 



R b /R = 0.710 - 0.716 (6) 
c b = (0.221 - 0.225) Mm/s 



The indicated range for Kb has been confirmed recently |22] 



Actually, within the present uncertainties, the information on R& and on q, are not 
independent, since to the per cent level, Cb and R& are related through 
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2GMg /Rq _ 1 j ^ 

3 Rq V Rft / 

Several determinations of the helium photospheric abundances have been derived from 
inversion (deconvolution) of helioseismological data, yielding (see Table |TV| ) : 

Y photo = 0.233 - 0.268 . (8) 

One should note that the small errors often quoted generally reflect the observational fre- 
quency errors only. The results actually depend on the method of inversion and on the 
starting physical inputs, e.g. the EOS. Such a dependence has been recently studied in 
Ref. ||. 

With these additional constraints, Eqs. (^) and (^), there are essentially no 
free parameters for SSM builders. 



C. Results of standard solar models 

Table | shows the main results of recent solar models produced by various authors; all of 
them include microscopic diffusion of helium and heavy elements, using (slightly) 
different physical and chemical inputs, which are summarized in Table ||. 

By varying Y and a, and adjusting (Z/X) in within the observed range, all models are 
able to reproduce L , R and (Z/X) p / loto at the solar age. 

The comparison with helioseismological measurements, which provide additional con- 
straints, is interesting, see Fig. |3| from |55 |. Only some of the recent solar models with 



microscopic diffusion are in agreement with the helioseismological constraints. With the 
notation of Table | they are P94 pf, RVCD96 |2|, BP95 0] and FRANEC96 @. On 



the other hand, all models without diffusion fail. The importance of diffusion for achieving 
agreement with helioseismological data is discussed in appendix A. 

We shall consider as Standard Solar Models (SSMs) only those models which 
pass the helioseismological tests. We recall that, among these, we shall refer to BP95 



"best model with diffusion" [ |15|| as to the Reference Solar Model (RSM). 

From Table | one sees that calculated central temperatures are in agreement at the per 
cent level. 

As a matter of fact, this shows that building up a stellar model is by now a well- 
established and reliable procedure, the small differences being essentially the result of small 
variations in the assumed input physics. At the same time, this agreement indicates the 
consensus about the treatment of the rather sophisticated physics needed to account for the 
behaviour of stellar structures. 

Concerning neutrinos, the solar model determines the internal distribution of tempera- 
ture, density and H abundance and, thus, the efficiency of the various chains. In this way 
one gets the fluxes at earth of the different neutrino components ($ pp , $ P ep, ^Be^ • • •)• From 
Table [QJ one notes that: 

i) the fluxes of pp neutrinos are stable to a few percent. As they are the most abundant 
($pp pa $tot) by far, their flux is directly related to the solar luminosity. 

ii) The ratio 
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f = $pep/($pep + $ PP ) 



(9) 



is almost independent of solar models, as it is weakly sensitive to the solar temperature and 
density. 

iii) After the pp, the 7 Be neutrinos are the most abundant, their flux accounting for about 
8% of the total, and again the prediction is rather stable among the different calculations 
with the exception of DS96. 

iv) The spread of the calculated 8 B neutrino fluxes is much larger showing the larger 
sensitivity to the different physical inputs used in different solar models. 

v) Concerning neutrinos from the CN cycle, again their prediction is somehow model- 
dependentfj. 

For a full CN cycle one should clearly have the same number of 13 N and 15 neutrinos. 
Actually, the reaction 14 N (p,7) 15 is too slow in the solar region below 10 7 K and the chain 
has not reached equilibrium, still favoring the transformation of 12 C into 14 N . This is why 
13 N neutrinos are slightly more abundant than the 15 neutrinos, i.e., if one defines 

r] = $ N /($ N + $ ) , (10) 



one finds r\ slightly larger than 0.5. Note that r) = 0.53 in most calculations. Table |TJ does 
not include results on 17 F neutrinos, which are two orders of magnitude less abundant than 
13 N or 15 O neutrinos^. 

We find useful to group some neutrinos according to their origin, putting pp and pep 
neutrinos in the same group 

$ P = $ PP + $ P ep , (11) 
and into another the neutrinos from the CN cycle: 

$CNO = $N + $0 • (12) 

Following the common terminology, we call these latter the "CNO neutrinos" , although they 
are neutrinos just from the CN cycle. 

Figure |5] shows the production region in the sun for the main neutrino components. 
Because of the strong temperature dependence, 8 B neutrino production is peaked at a very 
small distance from the solar center (R^0.04 R©). The same holds for 13 N and 15 neutrinos, 



3 In |32| Shaviv presented values of <1?n anc l an or< ler of magnitude smaller than those of all 
other calculations. This feature is no more present in the more recent calculation by the same 



author [23] 



4 The reader will note that the model DS96 ]^3| is the one yielding the smallest prediction on 
7 Be , CNO and 8 B neutrino fluxes (Table [Hl|) and the smallest central temperature. We remark 
that this model is not a standard solar model according to our definition, since it does not satisfy 
the helioseismological constraints. In addition the astrophysical S-factors used in [23] seem to us 
not acceptable, see Sects. [V B| and [V|. 
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which are produced with equal rate up to Rps0.14 R q . The secondary peak at R^0.16 R 
is due to 13 N neutrinos only, in a region where the CN cycle is marginally active, C and N 
have not reached their equilibrium values and nuclear reactions are transforming the more 
abundant 12 C into 14 N. The 7 Be neutrino production peaks at R~0.06 R Q , whereas for 
pp+pep neutrinos it peaks at R^O.l R Q and extends over a large portion of the solar core. 



D. Solar neutrino spectrum and the predicted neutrino signals 

The energy spectrum of each component (pp, pep, 7 Be , 8 B , . . . ) is deter- 
mined by kinematics and/or nuclear physics and it is (essentially) independent 
of solar physics. 

The pp neutrinos have a continuous energy spectrum extending up to E v = 0.42 MeV 
||13|| . For a pep neutrino, its momentum has to match that of the produced deuteron and 
consequently its energy is uniquely determined by kinematics, giving E v = 1AA2 MeV. The 
energy of 7 Be neutrino depends on the state of 7 Li produced in the electron capture reaction. 
One has E v = 0.861 MeV with 90% probability ( 7 Li ground state), and E u = 0.383 MeV 
with 10% probability (first excited state). The energy spectrum of 8 B neutrinos extends up 
to E v = 15 MeV |I3| , |57[ . Apart for the very rare hep neutrinos, 8 B neutrinos are thus the 
most energetic ones. In the CN cycle, the neutrinos from (3 + decay of 13 N and 15 have end 
points at 1.20 and 1.73 MeV, respectively. 

In Table [V] we give the average energy for each neutrino component. The average energies 
for the p and CNO neutrinos have been calculated by taking the values of rj and £ as in the 
RSM. 

The solar neutrino spectrum predicted by standard solar models is shown in Fig. |] 
(from Ref. ||24|| ). Note that the pp neutrinos, which are the most abundant and the most 



reliably predicted ones, are also those with the lowest energy and thus are most difficult to 
detect, since neutrino cross sections increase with neutrino energy. On the other hand, 8 B 
neutrinos, which are rare and delicate to predict, are those with the highest energy so that 
their detection is comparatively "easy" . 

The solar neutrino signal in radiochemical experiments is expressed as the reaction rate 
per target atom. A suitable unit is the Solar Neutrino Unit (SNU), defined as one reaction 
per second per 10 36 atoms. A 1 SNU signal means that in a target containing 10 31 atoms (of 
the order of several hundreds tons) there is about one reaction per day; this already gives an 
idea of the low counting rate of these radiochemical experiments. The signal S is expressed 
in terms of the neutrino interaction cross section a(E u ) and the differential neutrino flux 
d<&/dE v : 



S= / dE v a{E v )d$/dE v , (13) 

JE t h 

where E t h is the energy threshold for the detection. In the case of standard neutrinos, the 
spectrum can be written for each component as 

d$ i /dE v = $ i df/dE v , (14) 
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where df/dE u is normalized to unity and it is fully determined from nuclear physics and 
kinematics. The signal can thus be expressed as the sum of contributions arising from each 
neutrino component 

S = J2°i®i, (15) 

i 

where 

<Ti = J dE u a(E u )df/dE u (16) 

is the detection cross section averaged over the spectrum of the i-th neutrino component, 
and is again determined purely from nuclear physics. 

In summary, nuclear physics is responsible for cr,j, whereas solar physics de- 
termines $j. 

Of course cr, depends on the reaction which one is considering and we present in Table [V] 
the quantities relevant for the two reactions so far used in the experiments 

u e + 37 CI -> e~ + 37 Ar ( E th =0.8U MeV) (17) 

and 

u e + 71 Ga -> e~ + 71 Ge (E^=0.233 MeV) , (18) 

where the threshold energies are indicated in parenthesis^]. 

The cross sections increase with the neutrino energy. This is a general feature of weak 
interaction cross sections at energies well below the Fermi energy scale. In addition, this 
trend is more evident at low energies due to the presence of reaction thresholds, so that 
only neutrinos with high enough energy can induce the process. Note that the 8 B neutrinos, 
which are the most energetic ones, have cross sections three order of magnitude larger than 
the others, and this compensates for their significantly smaller flux. 

The contributions (o"j$j) to the signal in the Chlorine and Gallium experiments from 
each component of the neutrino flux are shown in Table |V1J], as predicted by the RSM. For 



Chlorine most of the signal comes from 8 B neutrinos, whereas for Gallium the largest contri- 
bution is from pp neutrinos. This is a consequence of the quite different energy thresholds 
for the two reactions. 



Table |T| shows the total signals in the Chlorine and Gallium detectors, 

S G a = J2 a Ga,i$i Sci = J2 °Cl,i®i , (19) 

i i 

as predicted by the different SSMs. 



5 Note that and v T cannot induce reactions ([l?]) and ( |l8| ) as the corresponding charged lepton 
{ji or r) would require neutrino energies orders of magnitude larger than those appropriate for the 
sun. 
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For the Chlorine signal the predicted values vary among the SSM within 20%, mainly 
due to the variations of the 8 B flux. The spread of predictions for the Gallium experiments 
is smaller, reflecting the rather stable estimates on pp and 7 Be neutrino fluxes. 

Radiochemical experiments are not the only way to detect neutrinos. As an alternative 
approach one can use the neutrino-electron elastic scattering 

v + e~^v + e~ , (20) 

detecting the scattered electron. 

Whereas for neutrino interactions with nuclei the cross section depends on nuclear physics 
and it is thus subject to some uncertainties, the cross section for scattering on electrons is 
well known from elementary particle physics. Both charged and neutral current interactions 
occur for electron neutrinos. For and v T only the neutral current interaction contributes 
and this results in a reduced cross section, e.g. for KAMIOKANDE o v ~ (l/7)a Ve . If one 
detects electrons with energy E e larger than some minimal value E Q , the signal is: 

S(E e >E o) =f iEAE ^f^*}^ , (21) 

where da/dE e is the differential neutrino-electron elastic cross section, which depends on 
the energy of the scattered electron E e . Because of background limitations, E Q is generally 
larger than a few MeV, so that only 8 B neutrinos are detected. 

Assuming standard neutrinos, the signal can thus be translated into an ef- 
fective (energy integrated) 8 B neutrino flux, and data are generally presented in that 
form. Note, however, that this "experimental flux" can be taken as the true 8 B neutrino 
flux only for standard neutrinos: should v e transform into u^, the detection cross section 
would be different. 



E. A look at experimental results 



So far we have results from four solar neutrino experiments; see Table VI for a summary 



and Refs. p^ , p5| for detailed reviews. 

The KAMIOKANDE experiment || , located in the Japanese Alps, detects the Cerenkov 
light emitted by electrons that are scattered in the forward direction by solar neutrinos, 
through the reaction (POP- 

This experiment, being sensitive to the neutrino direction, is the prototype of neutrino 
telescopes and is the only real-time experiment so far. 

The experiment is only sensitive to the high energy neutrinos [E v > 7 Mev) from 8 B de- 
cay. The solar neutrino spectrum deduced from KAMIOKANDE is in agreement (within 
large uncertainties) with that of neutrinos from 8 B decay in the laboratory ||. Assuming 
that the spectra are the same (i.e. standard z/ e ), one gets for the 8 B neutrino flux the result 
shown in Table [VI|. This corresponds to about 500 neutrino events collected in a 5.4 year 
data-taking period. Although there is no really sound way of combining statistical and sys- 
tematical errors, we use the usual recipe and combine them quadratically (see next section); 
here and in the following , we take 
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$ XA = (2.73 ± 0.38) • 10 6 cm- V 1 



(22) 



This value is less than one half that predicted by the RSM. 

All other experiments use radiochemical techniques. The 37 CI experiment of Davis and 
coll. H has been the first operating solar neutrino detector. The reaction used for neutrino 
detection is the one proposed by Pontecorvo in 1946 ||1|: 

v e + 37 CI -f e~ + 37 Ar . (23) 

The energy threshold being 0.814 MeV, the experiment is sensitive mainly to 8 B neutri- 
nos, but also to 7 Be neutrinos. The target, containing 10 5 gallons of perchloroethylene, is 
located in the Homestake gold mine in South Dakota. Every few months a small sample of 
37 Ar (typically some fifteen atoms!) is extracted from the tank and these radioactive atoms 
are counted in low background proportional counters. The result, averaged over more than 
20 years of operation, is 0] 

Sci = 2.55 ± 0.25 SNU. (24) 

The theoretical expectation is higher by a factor three. For almost 20 years this discrepancy 
has been known as the "Solar Neutrino Problem" . About 75% of the total theoretical 
rate is due to 8 B neutrinos, see Table |V1J| , and hence it was for a long time believed that 
the discrepancy was due to the difficulty in predicting this rare source. 

Two radiochemical solar neutrino experiments using 71 Ga are operating: GALLEX, lo- 
cated at the Gran Sasso laboratory in Italy and using 30 tons of Gallium in an aqueous 
solution, and SAGE, in the Baksan valley in Russia, which uses 60 tons of gallium metal. 
The neutrino absorption reaction is 

v e + 71 Ga -> e~ + 71 Ge . (25) 

The energy threshold is E = 0.233 MeV, and most of the signal arises from pp neutrinos 
with a significant contribution from 7 Be neutrinos as well (see Table |VU|) . In each experi- 
ment, the rate of neutrino interactions in the Gallium tank is about 0.5 event per day. The 
Germanium atoms are removed chemically from the Gallium and the radioactive decays of 
71 Ge (half-life=11.4 days) are detected in small proportional counters. The results of the 
two experiments, see Table [VT|, can be combined (see next section) to give 

S Ga = 74 ± 8 SNU, (26) 

which we shall use as the representative value of the Gallium signal. Again the value is 
almost a factor two below the theoretical prediction. 

An overall efficiency test of the GALLEX detector has been performed by using an intense 
51 Cr neutrino source, (61.9 ± 1.2) PBq 00. The source, produced via neutron irradiation of 



6 The preliminary results on 51 Cr neutrino source measurements by SAGE collaboration have been 
presented recently in ||. 
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about 36 Kg of cromium enriched in 50 Cr, primarily emits 756 keV neutrinos. It was placed 
for a period of a few months inside the GALLEX tank, to expose the gallium chloride target 
to a known neutrino flux. The number of observed neutrino events agrees with expectation 
to the 10% level. This result "(a) provides an overall check of GALLEX, indicating that 
there are no significant experimental artifacts or unknown errors at the 10% level that are 
comparable to the 40% deficit of observed solar neutrino signal (b) directly demonstrates 
for the first time, using a man-made neutrino source the validity of the basic principles of 
radiochemical methods used to detect rare events ..." 0. And last but not least: "because 
of the close similarity in neutrino energy spectra from 51 Cr and from the solar 7 Be branch, 
this source experiment also shows that the Gallium detector is sensitive to 7 Be neutrinos 
with full efficiency.' 1 fTj. 

F. A remark on errors 

Because of the subtleties of experiments, of the difficulty in estimating errors of theoret- 
ical evaluations and in face of the persistent discrepancy between experimental results and 
standard theory (standard solar models and standard neutrinos), the discussion about the 
meaning of errors is lively at almost any conference on solar neutrinos. In this section we 
intend to define the errors we adopt, describe what we want to do with these errors, provide 
some justification to our choice and outline its limits. 

Concerning experimental errors, we add in quadrature statistical and systematical errors. 
We tentatively use this global error as an indicator of the distance between the true and the 
measured values of a physical quantity, according to the rules of Gaussian statistics. 

Essentially we are assuming that the systematical uncertainty is a variable which can 
fluctuate according to statistical rules, so that we consider statistical and systematical errors 
as independent of each other and on the same footing. This choice seems to be crazy (and 
actually it would be in many cases), however it has some justification in the present context. 

As an example, consider measuring the period of a pendulum using a watch. This watch 
is guaranteed to measure time better than one second per day: one does not know if it 
is fast or slow, but it does not gain or loose more than ls/day. This systematical error 
cannot by traced to a single origin (otherwise it would be eliminated), but arises from many 
different sources, essentially the mechanical tolerances on each of the many components of 
the watch. Each of these sources gives a small contribution, positive as well as negative, to 
the final error. Also each contribution is statistically independent from the others. Under 
these conditions one can reasonably suppose that, if one produces a large number of watches 
of the same kind, the systematical errors are distributed according to a gaussian with null 
average, i.e. the distribution is the same as for statistical fluctuations. 

Clearly one does not gain any accuracy by repeating the measurements over several 
oscillations, in exactly the same conditions with the same watch. On the other hand if 
one is using N different watches of the same kind (each one with its systematical error 
e Q = ± ls/day) and combine the results, accuracy is improved, as the chance of having a 
fast and a slow watch are equal. In other words, in this example systematical errors can be 
treated similarly to statistical fluctuations. 

The situation, at least in radiochemical experiments, is somehow similar. For instance 
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the 33 runs of GALLEX are actually 33 independent experiments: after each run, the 
71 Ge atoms are extracted in specific conditions, which are independent of the run number, 
and placed in a new counter. The radon content of the atmosphere where extraction occurs 
fluctuates around some mean value, within some systematical uncertainty. All counters have 
approximately the same efficiency, again within some systematical uncertainty, but one is 
better and another is worse than the estimated mean efficiency Thus the systematical 
error of the full experiment is significantly reduced with respect to that pertaining to each 
individual run and, it can approximately be determined by statistical considerations. 

The same procedure used to combine different runs can be used in principle to combine 
two different experiments like SAGE and GALLEX. The resulting systematical error is 

l/e 2 {syst) = l/e GALLEX (syst) + l/e 2 SAGE (syst) , (27) 
in analogy to the statistical one 

l/e 2 (stat) = l/e 2 GALLEX (stat) + l/e 2 SAGE {stat) ; (28) 
and the global error on Gallium result becomes 

e 2 = e 2 (syst) + e 2 (stat) . (29) 



This is how Eq. (^6|) has been obtained. 

The situation is more complex when considering errors of theoretical results. When these 
errors correspond to uncertainties on some measured quantities which are used as input in 
the theoretical calculations, we are back to the previous case. Often, however, additional 
uncertainty is added due to extrapolations of experimental data, e.g. when estimating nuclear 
cross sections at energies of astrophysical interest. In other cases (e.g. opacity tables) the 
only information comes from theoretical calculations. In these cases one can only resort to a 
comparison among different calculations. A total theoretical range €ttr is defined by Bahcall 
as "the range in values of published state-of-the-art calculations" [E1,0], and this might be 



considered as the equivalent of an experimental three sigma error, i.e. the probability of 
being outside this range is about 0.3%. This looks plausible, how many sigmas being clearly 
a matter of taste. We generally adhere to Bahcall's choice. For consistency of notations, as 
we always show experimental la errors, we shall add to the theoretical quantities an error 

e theo = -tTTR ■ (30) 

All this clearly means that treating errors given by equations of the form ( p9|) 
or (|30|) as if they were indicators of statistical (Gaussian) fluctuations is just a 
rough approximation, or even an undue simplification. 

The confidence levels obtained in this way are thus to be taken cum grano salis, or much 
more. 
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II. NEUTRINO FLUXES ALMOST INDEPENDENTLY OF SOLAR MODELS 



All experiments give signals significantly smaller than those predicted by the RSM and 
by other SSMs, for standard neutrinos. Is this a problem of SSMs or is there anything 
deeper? Insight on this question can be gained by looking at properties of neutrino fluxes 
which are largely independent of solar models. In essence, one is asking if experimental data 
are in agreement with the following assumptions: 

A) the solar luminosity is supported by H burning reactions; 

B) electron neutrinos do not disappear in their travel from sun to earth (i.e. standard 
neutrinos) . 

The first assumption looks rather innocent, the second one being really the hypothesis to 
be tested. We shall follow this approach by discussing: a lower bound on the Gallium signal 
for standard neutrinos (Sec. [II B|) ; the relationship between KAMIOKANDE and Chlorine 
signals (Sec. [11 C| ); the information on intermediate energy neutrinos derived by using the 
full set (or any subset) of the data (Sec. [II D|) . 



We shall show that the experimental results look mutually inconsistent if the above 
two assumptions hold true. The point is that the flux of intermediate energy neutrinos 
(E u » 0.5-2 MeV): 

$int = $ Bc + $CNO + $ P ep (31) 

would have to be negative! 

We shall see that the same occurs if we disregard any of the experiments. Actually we 
do not know of any sound argument to doubt the experimental results, but we 
shall entertain this possibility, just as a working hypothesis, to test the consistency 
of the assertion: neutrinos are standard and (some) experiments are correct. We would like 
to make this point as clear as possible: when saying "disregard one experiment" or "assume 
it is wrong" we are advocating some hypothetical, large, unknown systematic error, just for 
the sake of saving/testing the hypothesis of standard neutrinos. This is clearly a desperate 
hypothesis, at least for 51 Cr-calibrated Gallium detectors. 

The inconsistencies will be derived by treating the different components of the neutrino 
flux as essentially free parameters. Additional information can be gained from the physical 
relationships among the neutrino components, as imposed by the development of the fusion 
chains. For example, pep neutrinos accompany pp neutrinos, with a proportion which has 
to be essentially the same in any solar model. Also, if the CN cycle is efficient, there 
will be about as many 13 N as 15 O neutrinos. By using these additional, however weak, 
assumptions one gets more severe hints against the hypothesis of standard neutrinos and/or 
tighter constraints on 7 Be and CNO neutrino fluxes (see Sees. |II E| and [II F| ). Quantitative 



statements about the chances of standard neutrinos and upper bounds on intermediate 
energy neutrino fluxes are presented in Sec. [11 G| . These bounds appear in conflict with the 
results of SSM, particularly for 7 Be neutrino flux (Sec. |IIH|) . 



Futhermore, since 7 Be and 8 B neutrinos both originate from 7 Be nuclei, there is a rela- 
tionship between them which looks again in contradiction with the experimental data, for 
standard neutrinos (Sec. [1T1| ). 

All these (essentially) solar model independent arguments point towards some non- 
standard neutrino properties, the main ingredient being the luminosity constraint A). In 
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Sec. H J| we relax even this hypothesis, still finding indications towards non-standard neutri- 
nos. On the same grounds, we discuss in Sec. [II K| the case of universal neutrino oscillations 
in the framework of a solar model independent approach. 

Before presenting the discussion let us advance a warning. We regard this section as the 
most relevant one of the present review and we therefore tried to be as detailed as possible. 



For a first reading we recommend Sees. II A and II B, where the basic ideas are explained, 



and Sees. [11 h] and |11 FJ , which contain the main results, and the concluding remarks at the 



end of this section ||. 



A. The solar luminosity constraint 



The solar luminosity constraint, essentially given in Eq. (0), will be a major ingredient 
of this section, so that it is important to express it more precisely and to discuss it accuracy. 

Assuming that all solar energy originates from nuclear reactions, and that all these reach 
completion, i.e. : 

4p + 2e~ -> 4 He + 2u e , (32) 
one immediately gets the following constraint for the neutrino fluxex $j: 

K Q = £ (f - (E v )i) ^ , (33) 

where Q = 26.73 MeV and (E v )i is the average energy carried by z-th neutrinos. The 
coefficients 

0* = f- (Euh (34) 

represent the average electromagnetic energy released per emitted neutrinos, and can be 
calculated from Table [V|. Except for the rare 8 B neutrinos, (E u )i is of the order of half MeV, 
so that Eq. (g) — corresponding to Eq. (^) when {E v )i is neglected — is accurate to the 
five per cent level. 

Concerning the accuracy of Eq. (|33|), the following comments are in order, see [|58|j ]: 



i) gravitational energy generation is neglected. According to standard solar model cal- 
culations, this causes an error of about 3 ■ 10~ 4 | J13[| , negligible in comparison with the 
uncertainty on K®. 

ii) The abundance of 3 He nuclei is assumed to be in equilibrium, which is not strictly 
correct. In the outer region of the solar core 3 He is continually produced, but the temperature 
is too low to burn 3 He at the equilibrium rate. Thus the pp chain is not always completed. 
Similar considerations hold for the CN cycle. All this gives, according again to standard 
solar model calculations, corrections to Eq. (^) of the order of 4-10 -4 [58 . 

All in all, Eq. (|3"3"P should be accurate to better than 1%. 



7 In the same paper one derives direct upper bounds on neutrino fluxes by using Eq. (33) 
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Finally, we remark that one assumes the sun to be at thermal equilibrium when using 
Eq. (33) to relate the present solar constant to the present nuclear energy production 
rate and present neutrino fluxes. The assumption of thermal equilibrium (stationary sun) 
looks quite reasonable to us. For the reader who is willing to abandon it, we have anyhow 
prepared Sect. [11 J| . 



B. The minimal Gallium signal for standard neutrinos 

For standard neutrinos, one can find a lower bound to the Gallium signal by the following 
considerations |59|j : 

i) Since the pp neutrinos are the least energetic ones, they have the smallest cross section 
in Gallium detectors, i.e. <JGa,i > Ga,pp- Thus from Eq. ( |T9"D one derives 

Sea > CrGa,pp^2®i = ^Ga,pp^tot , (35) 

i 

i.e. the minimum signal is obtained by assuming that all neutrinos are from the pp reaction. 

ii) The solar luminosity, namely the electromagnetic energy released per unit time, is 
obtained with a minimum number of fusions when neutrinos carry away the least energy. 
The minimum total flux $ m j n is thus obtained when all neutrinos are from pp, and it is thus 
related to the solar constant K & : 

$tot > $min = , (36) 

~B pp 

where the term in the denominator is the average electromagnetic energy released per emitted 
pp neutrino. This gives thus: 

K 

^Ga — °Ga,ppJj- ■ {<$ < ) 

^%pp 

In other words, the minimum Gallium signal is obtained when all neutrinos are from the pp- 
reaction. Any solar model will give a prediction not smaller than Sg* n , provided 
only that the present electromagnetic energy production rate in the sun equals 
the 'presently observed solar luminosity. 

By using the values previously given K Q = (0.853 ± 0.003) • 10 12 MeV cm -2 s _1 and 
(?Ga,pp = (1.18 ± 0.02) • 10- 9 SNU cm 2 s , one derives: 

gmm = 77±2 g NU; ( 3g ) 

where the error essentially comes from the uncertainty on <JGa,pp- The central value of this 
prediction is already 3 SNU above that of the experimental result (74±8 SNU), however the 
distance between the two values is well within the errors. Gallium results do not violate the 
lower limit for solar neutrinos. 

However, this is far from being a satisfactory result. If pp burning is at work, it appears 
difficult to switch off the pep reaction. The pep neutrinos being more energetic, the minimal 



signal clearly increases. In Ref. the value Sg* n = 79 SNU is obtained under the physically 
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very plausible assumption that together with pp neutrinos also pep neutrinos are created, 
with the ratio £ = $ pe p/ ($pp + $pep) ~ 2.5 • 1CT 3 as given by the solar model of the same 
Ref. ||59|| . We did not consider pep contribution as we wanted to use minimal assumptions; 
should we follow the approach of Ref. |59|], with the cross sections given in Table [V| and 



£ = 2.36 ■ 1(T 3 we would find, in place of Eq. S^ n = 80 ± 2 SNUQ, which is still 

consistent with experimental result at la level. 

However the situation gets somehow puzzling if one takes into account the other mea- 
surements. For example, KAMIOKANDE does observe the energetic 8 B neutrinos, implying 
a 8 B contribution to the Gallium signal of about 7± 2 SNU, for standard neutrinos. All this 
means that for standard neutrinos the signal in Gallium experiments has to exceed 87±3 
SNU. The difference between this minimal expectation and the measured Gallium signal 
gets now larger than the error. 

By considering the above information, one finds that the contribution of 7 Be and CNO 
neutrinos is: 

S B e+CNO = (-13 ±8.5) SNU. (39) 

Therefore, it should not exceed 13 SNU, if the observed signal has to agree with with the 
minimal estimate to the three sigma level. This is significantly smaller than the predictions 
of all SSMs (about 50 SNU, including 38 SNU arising from 7 Be neutrinos, see Table |VU| ). 

In conclusion, Gallium results are still consistent with standard neutrinos, essentially 
implying that the vast majority of them are from the pp reaction and that intermediate 
energy neutrinos are definitely fewer than in the SSMs. 

The recently proposed Gallium neutrino observatory GNO KJ which foresees a hundred 
ton target and improved detection techniques, aims at a measurement of the Gallium signal 
with an accuracy of 5 %. This would provide a significant progress in the "minimal Gallium 
signal" test. 



C. Chlorine and KAMIOKANDE 

As soon as the KAMIOKANDE data were available, assuming standard neutrinos, it 
became clear that they presented a new puzzle, when compared with the Chlorine results. 
Let us divide the Chlorine signal in two parts, corresponding respectively to intermediate 
energy neutrinos and to 8 B neutrinos: 

Sc7 = Sci,int + $Cl,B ■ (40) 

KAMIOKANDE only detects high energy events (E v > 7 MeV); assuming standard 
neutrinos this is enough to tell the energy integrated flux of 8 B neutrinos, and thus their 
expected contribution to the Chlorine signal, Sqi,b, can be calculated, by using the estimated 
value of <tci,b (see Table |V|): 



The small difference with respect to the result of Ref. [59] originates from the slightly different, 
more recent values for OQa^-p-, &Ga,pep and £ which we are using. 
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Sci,b = vci,B<5>B am = 2.98 ± 0.41 SNU , (41) 

where most of the error comes from uncertainties on the KAMIOKANDE signal. By compar- 
ison with the measured Chlorine signal Sci = 2.55±0.25 SNU, one sees that the 8 B neutrinos, 
as measured by KAMIOKANDE, should produce a signal in Chlorine larger than observed. 
So little space is left for the intermediate energy neutrinos: 



~>CLint 



-0.48 ±0.49 SNU. (42) 



Note that the central value is negative, however Eq. fl42D is clearly consistent with zero. 
At "three sigma" level one has Scimt < 1 SNU, this limit being a factor two smaller than 



the RSM prediction, see Table VII 



Let us remark that the above estimate on Sci,int is independent of the stationary sun 
hypothesis A). 

The present uncertainty on Sci t mt depends essentially on the systematical error of the 
KAMIOKANDE result. A better determination of Sci,i n t would require a strong reduction 
of the systematical error in SUPERKAMIOKANDE with respect to KAMIOKANDE. In 
this case the ultimate uncertainty would be that due to the Chlorine experiment, implying 
ASci,mt ~ 0.25 SNU, which would provide significant information on intermediate energy 
neutrinos. 

D. The space for intermediate energy neutrinos 



An extension of the minimal Gallium signal argument of Sec. |TB| can elucidate the role 
of various experiments in providing solar model independent information on neutrino fluxes. 
We already known that experimental results are smaller than the expectation from SSM and 
consequently the fluxes of intermediate and high energy neutrinos should be smaller than 
theoretically predicted. The relevant question is: how large neutrino fluxes are allowed? 

We have three equations relating the solar constant, the Gallium and Chlorine signals 
with the fluxes $>f. 

K Q = J2Qi®i, (43) 

i 

S G a = Y. a Ga,i^i, (44) 

i 

S CT = °ci,i$i , (45) 

i 

and the 8 B neutrino flux as measured by KAMIOKANDE 

^Kam = ( 2 73 ± 0>38 ) . 10 6 cm -2 s -l ( 46 ) 

In order to determine the maximal allowed fluxes of intermediate energy neutrinos let 
us assume that only one flux $^ (k=Be, N, O or pep) is non vanishing and attribute to 
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$k all the signal in the Chlorine and GALLEX experiments pertaining to the intermediate 
energy neutrinos. 

For the Chlorine experiment, one has just to subtract from the experimental value the 
8 B contribution, i.e. the maximal flux $^ satisfies 



■>ci 



^CZ,fc$k + ^Cl,B^B • ( 47 ) 



For handling the information provided by Gallium experiments, we derive $ pp from the 
luminosity equation. In terms of the average electromagnetic energy Qi released per emitted 
neutrino, see Eq. B3, one has from Eq. 



*pp = ^-pr*B-E*# («) 

Vpp Vpp int ^PP 

so that 

Sea = T^O"Ga,pp + {<TGa,B ~ J^-^Ga,pp)^~Q + ^{^Ga,i ~ ^-^Ga,pp)^i ■ (49) 

UJpp Wpp int UJpp 

All terms in the brackets are of course positive as the minimum signal is obtained when 
all neutrinos are from pp. 

Again isolating the k-th component, one obtains the maximal allowed flux from: 

SGa = VGa,ppJ^- + ((TGa,B - T^- a Ga,pp)®B + ( a Ga,k ~ ^^ a Ga,pp)^k ■ (50) 
^ipp ^ipp ^ipp 

The reader recognizes, for the limiting case = $b = 0, the inequality found in the 
previous section. 

Clearly this exercise is most interesting for 7 Be neutrinos, as they are predicted to be 
the second most abundant ones. By using the cross sections in Table [V| and taking k=Be, 
from Eqs. ( f£T| , [50"D one gets: 

S OJ = 0.24$ Be + 1.11$ B (51) 

and 

S G a = 77 + 6.19$ Bc + 2.43$ B , (52) 

where the signals are in SNU, $g e in 10 9 cm -2 s -1 and $g in 10 6 cm~ 2 s _1 0. 

In this way the experimental information can be presented in the (^Bj^Bc) plane, see 
Fig. H Dashed lines are obtained by adopting in Eqs. (|5T|) and ([52]) the central values 
of experimental results. Full lines correspond to the experimental results ±lcr. Even by 
allowing for such a spread, still the curves intersect in the unphysical region ($Be — 0)- 



We note that errors on the cross sections are negligible with respect to the uncertainties on 
experimental signals. 
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In Fig. [|we also introduce the information from KAMIOKANDE experiment, Eq. (|46"D, 
and the following conclusions can be drawn: 

i) the la allowed areas intersect in the unphysical ($ Be < 0) region. The same holds if 
one disregards either the Chlorine or the KAMIOKANDE experiment. Should one disregard 
Gallium result, still at la one has a very small value $ Bc < 1.0-10 9 cm -2 s -1 . 

ii) in order to stay within 2a from each experimental result the beryllium flux has to be 
$ Bc < 1.5-10 9 cm -2 s -1 , i.e. a factor three at least smaller than the predictions from SSMs. 



Similar considerations hold for the other intermediate energy neutrinos, see Table VIII 



In any case, the best fit point is in the unphysical region and no point can be found in the 
physical region which is within la from each experiment. 

Roughly speaking, the hypothesis of standard neutrinos (which requires positive fluxes) 
disagrees with each of the three independent experimental results by at least la. The 
probability for such a situation is at most (0.32) 3 , i.e. a few percent. Little space seems to 
be left for intermediate energy neutrinos. 



E. Four equations and four unknowns 



So far the only assumption about the sun was Eq. fl43|) , that connecting neutrino fluxes 
to the solar constant. Stricter constraints on the intermediate energy neutrinos can be 
obtained by using some additional, albeit very weak, hypothesis. Theoretically, we know 
that pep neutrinos must accompany pp neutrinos and also that the CN cycle is (almost) at 
equilibrium when it is efficient in the sun. The precise values of £ = $pep/($pp + $pep) 
and rj = $n/($n + $q) are n °t important for the following discussion and we shall use 



f = 2.36 ■ 10~ 3 and r] = 0.53, according to the RSM [lj 



In this way we reduce our unknowns to just four variables: 

$p = $ PP + $pcp , $CNO = $N + $0 > $Be and $ B , (53) 

i. e. as many as the available pieces of information. With the numerical values in Table [V], 
from Eq. (f43|) (after dividing both sides by Q/2) and Eqs. ((44]) and (f4~5|) one gets: 

63.85 = 0.980$ p + 0.939<l> Be + 0.937$ C NO + 0-498 ■ 10~ 3 $ B (54a) 
S Ga = 1.23$ p + 7.32$ Be + 8.72$ C NO + 2-43$ B (54b) 
Sci = 0.38 • 10- 2 $ p + 0.24$ Be + 0.41$ CNO + 1.11$ B (54c) 

<5>>< am = $ B , (54d) 

where all fluxes are in units of 10 9 cm -2 s^ 1 but for the 8 B flux which is in units of 10 6 cm -2 
s _1 and signals of radiochemical experiments are in SNU. Uncertainties on the numerical 
coefficients (arising from errors on the estimated neutrino cross sections) have been omitted, 
since errors on the experimental signals are dominant. 
We thus determine the four fluxes: 
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$ p = 1.905a - 0.229S Ga - 1.56$f am + 83.1 
$ Be = -10.65 O{ + 0.57S Ga + 10.4$g am - 43.1 
$ C no = 8.6S cl - 0.33S Ga - 8.8$g am + 24.4 

$ B = . (55) 



inserting the signals given in Table |VI], one finds^j: 



$ p = (66.7±2.0)-10 9 cm- 2 s" 1 
$ Be = (0-4 ± 6.6)-10 9 cm' 2 s" 1 
$ CN0 = (-2. ± 4.8)-10 9 cm" 2 s- 1 . (56) 

It is remarkable that, for standard neutrinos and with minimal and quite reasonable 
assumptions about solar models (i.e. the values of £ and 77), the main components of the 
solar neutrino flux are fully determined from available experiments. In other words, a 
solar neutrino spectroscopy is already at hand and it could be used to study the 
solar interior, if we know enough about neutrinos. 

We see, for example, from Eq. (|56"D that the pp+pep flux is determined for standard 
neutrinos with an accuracy of about 3%. This result might be surprising when considering 
that all experiments have no more than 10% accuracy. The point is that the total flux is 
fixed by the luminosity, whereas the cross section depends crucially on neutrino energy, so 
that approximately (a Ga ^ Be » a Ga ^ p ): 

( A< M ~ (^Sca) °Ga, PP , 7) 

{ $p J ~ I S Ga J a Ga>Be ' 1 ! 

The result on $CN0 clearly shows — again for standard neutrinos — that the CN cycle 
cannot be the main energy source in the sun, otherwise one should have: 

IK 

$ CNO ~ "7^ ~ 60 ■ 109 cm " 2 s" 1 , (58) 

in violent contradiction with Eq. (|56]). If the CN cycle were to dominate solar energy 
production, then the Gallium signal would be 

IK 

S G a « -^PGa,CNO « 550 SNU , (59) 

an order of magnitude larger than the actual result. 



10 When errors on the neutrino cross section are taken into account, the errors quoted in Eq. (56) 
are slightly enlarged, becoming ±2.2, ±7.2 and ±5.1, respectively. 
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F. Where are Be and CNO neutrinos? 



One notes that the central value of $cNO m Eq. ( |56| ) is negative i.e. unphysical. In view 
of the estimated errore, this does not seem to be a problem. However, there is a strong 
correlation between $cNO and ^Be; so that if $CN0 * s forced to be positive, then $g c 
becomes negative, and vice versa . 

In order to understand what is going on, and to clarify the role of each experimental 
result, let us once more reduce the number of unknowns. We start again from the basic 
equations ([43|)-((46]) and use the following tricks, similar to those used previously: 

(a) we group the neutrino fluxes as in Eq. (|53|), so that we are left with the four variables 

$p,$Be,$CNO and $ B - 

(b) Since (E u )q^q > (E u )-Q e , the corresponding cross section is larger than that of 

7 Be neutrinos. Thus the minimal CNO signal is obtained with the replacements 

(^>CNO ( E »)Be and a CNO -> CTBe ■ (60) 



This corresponds to dropping terms containing $CNO m equations ( |54a| - |54c!D and replacing 



$B C with $g e + $CNO' so that only the combination $g e + $CNO enters. We remark that 
such a substitution represents also a safe approach, since the theoretical value of ac a ,Be has 
essentially been verified (to the 10% level) by the GALLEX neutrino source experiment [7|, 
whereas only theoretical predictions exist for ac a ,CNO- 

One can eliminate <3> p by using the luminosity equation ( 54a|) , and each experiment 
provides a constraint on $Be+CNO and/or $g: 

S Ga = 80.1 + 6.14$ Bc+C NO + 2.43$ B 
Sci = 0.248 + 0.236$ Bc+ CNO + LH$B 

= $ B . (61) 

The result of each experiment can be plotted in the ($g, $Bc+CNo) plane, as shown in Fig. 
With respect to the situation of Sec. OD| , the intersection moves towards even more 



negative values of "^Be+CNO see Table |V11J| , and the allowed region in the physical part 



of the plane shrinks. Whichever experiment is discarded there is no point in the physical 
region within la from the remaining two results. The region at 2a from each experiment 
now allows only "^Be+CNO < 1- This bound is stronger than we found previously, the main 
reason being that the small, but non negligible, contribution of pep neutrinos 
to Chlorine and Gallium signals is included from the beginning. 

The bound on the sum clearly holds separately for 7 Be and CNO neutrinos. For these 
latter, however, we can get tighter constraints by putting $b c = in Eqs. ( |54aH54d]) , see 



again Table VIII. In this way one finds that in the area within two standard deviations from 



each experiment $CNO — 0-7- 



G. Probabilities, confidence levels and all that 

The intermediate energy neutrino fluxes, which we derive from the experimental signals 
assuming standard neutrinos, favor negative, i.e. unphysical values. Here we shall discuss the 
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statistical significance of this information. We try to estimate quantitatively the probability 
of compatibility with standard neutrinos and to determine the significance of upper bounds 
on neutrino fluxes. 

For pedagogical reasons, let us discuss first a kindred situation where only one variable 
is present, namely the determination of the neutrino mass from f3 decay experiments. 



1. Neutrino mass and (3 decay 



Measurements of the electron energy spectrum from tritium decay determine a parameter 
m 2 (apparent squared mass), which a priori is just a combination of experimental data and 



theoretical inputs. In principle, m app can be positive as well as negative. This parameter can 
be identified with the actual u e squared mass m 2 only if the interpretation of the experimental 
data (underlying theory, estimated energy resolution, . . . ) is correct, and in that case it 
should be positive. Even so (m 2 pp > 0) experimental results rh 2 app might come out negative 
due to statistical fluctuations of the physical quantities entering the definition of m 2 app . 
The latest edition of the Particle Data Book |)T| considers four tritium f3 decay experi- 



ments, all giving m 2 app < 0, and presents the weighted average: 

m 2 pp = (-54 ± 30) eV 2 . (62) 
This result prompts the following questions: 

• (a) Is the interpretation of the data ("working hypothesis") correct? 

• (b) Which upper bounds can be set on the actual neutrino mass? 



In the discussion we follow — to some extent — the approach of Ref. Bl . We assume 



that the experimental results have a Gaussian distribution, centered at an unknown value 

i 2 

app' 



m 2 , the width being given by the quoted experimental uncertainty A = 30 eV 2 . 



(a) If the theory is correct (m 2 pp = m 2 ) then of course m 2 app > 0. Evidently, the 
probability of obtaining a value equal or smaller than m 2 pp = —54 eV 2 is maximal 
when the neutrino mass is zero (m 2 = m 2 =0) and in that case it is: 

^°-° 36 - <63) 



For m 2 pp > the probability is even smaller. Furthermore the requirement m 2 pp > 
is just a necessary condition for identifying it with the physical neutrino squared mass. 
Thus our answer to question (a): there is at most a 3.6% probability that the 
working hypothesis is correct. 

It is somewhat controversial how to set bounds on an observable when the experimental 
results lie outside the physical region, and more than one answer to our question (b) is 
possible: 
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• (bl) The first approach we consider is similar in spirit to that just presented. We can 
determine an upper bound P2 on the probability that the result is equal to or smaller 
than rh 2 app if the apparent squared mass exceeds a value m 2 ,: 

f (rh 2 app -m 2 )/A , 3.2/2 

Conversely, by fixing a Confidence Level C.L.= 1 — P2, we can determine the corre- 
sponding value of m 2 , such that the above equation holds. For instance, if we choose 
P2 = 1%, we find that at the 99% C.L. the apparent neutrino squared mass is 
less than 15.8 eV 2 . 

The interpretation of this limit is the following: answering question (a) we used the fact 
that only positive values are physical. Now, if we have some additional information 
telling us that m 2 > 15.8 eV 2 , then there is at most a 1% chance that the working 
hypothesis is correct. 

Note that the identification m 2 app = m 2 is not necessary a priori and actually it is one of 
the hypothesis being tested. Also, to reach our conclusion we needed only the probability 
distribution function (p.d.f.), f(a\a), giving the probability of observing an experimental 
result a if the true value is a. 

We remark that / is different from the distribution defining the probability that, given 
an experimental result a, the true value is a. This latter p.d.f., which we denote by g(a\a), 
is the one we are going to use in a second approach to question (b), the so-called Bayesian 
approach to confidence limits, described extensively in Ref. |61| and references therein. 



(b2) Here one assumes that the working hypothesis is correct and one seeks information 
on the mass m from the knowledge of experimental result m 2 . To be precise, we list 
the assumptions: 

(i) the measured quantity can be identified with the physical neutrino mass, m 2 pp = m 2 . 

(ii) Bayes' theorem holds: 



g{m 2 \m 2 ) = /f^ ^ ' , (65) 

J j{m z \m z ) 7r(m z ) am z 

where 7r(m 2 ) is the a priori p.d.f of the neutrino squared mass, which we define as: 

n ( m 2) = { 1 if m2 > (66) 
I otherwise. 



As we have chosen / to be Gaussian, g(m 2 \fn 2 ) is also Gaussian, but its m 2 domain 
is restricted to the positive axis. For a confidence level 1 — P3, the upper limit m 2 , to 
the neutrino squared mass is now given by 

p > = ■ (67) 
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For instance, if we choose P3 = 1%, this time we find = 47.4 eV 2 and we can 
say that at the 99% confidence level the neutrino squared mass is less than 
47.4 eV 2 . 

In the present context, for the same m 2 , one has: 

P 2 = Pi-P 3 ■ (68) 

Thus P2 is always smaller than P3, and correspondingly the bound on the mass found with 
the second approach is always weaker than the first one. 

We may say in conclusion that the two approaches answer two different questions: 

• (bl) assuming that the physical mass squared is larger than some value m 2 ,, what is 
the chance that the interpretation of experiments is correct and that the result does 
not exceed —54 eV 2 ? 

• (b2) Assuming that the interpretation is correct, what is the probability that the 
physical mass squared is larger than m 2 , and that experimental result does not exceed 
-54 eV 2 ? 

Correspondingly, at a given confidence level, one associates two different bounds on the 
squared neutrino mass, which correspond to two different physical assumptions and atti- 
tudes: the first case (bl) is essentially a way of testing the interpretation of the experiment, 
if one is confident in a minimal value for neutrino mass; the second one (b2), is a way of 
testing a theoretical prediction on neutrino mass, if one believes in the interpretation of the 
experiment. 



2. Statistics and the solar neutrino problem 

The analogy with the previous case should be clear now. In the language of the previous 
section, the right-hand sides of Eqs. ([55]) define four apparent neutrino fluxes $p PP , , $g P , 
and ^qnO' wn i° n can be identified with the physical fluxes assuming standard neutrinos. 

Different attitudes are possible and correspondingly different questions can be raised. If 
one wants to test the chance of standard neutrinos with minimal assumptions about solar 
physics, then the relevant question is: 

• (a) what is the probability that all the apparent fluxes are non negative? 

If one is confident in some solar models, yielding definite predictions/lower bounds on the 
fluxes, then again as a test of standard neutrinos the question is: 

• (bl) what is the probability that neutrinos are standard if the true fluxes are at least 
$>L,i, i n face of the available results? 

On the other hand, one who believes in standard neutrinos and wants to test solar models 
will be interested in upper bounds $[/,i on the fluxes. His question is now: 
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• (b2) assuming standard neutrinos, what is the probability that the true fluxes do not 
exceed $u,i, m view of the experimental results? 

As compared to the (3 decay, the only complication is that we deal now with several 
variables (fluxes) instead of just one (the neutrino squared mass). 

For this reason we resorted to Monte Carlo techniques. We generated a large ensemble of 
sets of four simulated signals (Gallium, Chlorine, KAMIOKANDE and the solar luminosity); 
each simulated signal was extracted by a Gaussian distribution with mean and width equal 
to its actual experimental central value and error. When appropriate we only considered the 
subsets of three simulated signals where one of the neutrino experiments (Gallium, Chlorine 
and KAMIOKANDE) was in turn excluded. Given any set of four (or three) simulated 
signals out of the ensemble, we derived three apparent fluxes: <3> p , $b and either <3>B e or 
<&CNO assuming the other one ($CNO or ^Bc) equal to zero. We have chosen not to derive 
from the simulated signals both $g e and $CNO> since the possibility of separating the two 
signals is critically dependent on the difference between the ratio of the average cross sections 
for Be and CNO neutrinos in the Chlorine experiment and the same ratio in the Gallium 
experiments. Moreover, the limit obtained by assuming either $g c or $CN0 to be zero is 
more conservative. Note that when we assume $CN0 = 0, the upper limit on $b 6 is in 
fact an upper limit for $B e and for $Be+CNO being the cross section for CNO neutrinos 
larger than the one for Be neutrinos. In practice the system of equations giving the signals 
as functions of the fluxes was inverted with the techniques of singular value decomposition 
that automatically take care both of the case of almost degenerate equations (if one decides 
not to combine the two Gallium data) and of the case of an over determined system (when 
the apparent fluxes are less than the number of equations), giving in this last case effectively 
the best x 2 fit- In summary, we have generated an ensemble of a milion apparent fluxes 
from the ensemble of a milion simulated experiments. If $g e is positive in this ensemble 
60000 times out of a milion, we say that the probability of standard neutrinos is 6%. The 
other limits are derived similarly considering different experiments and/or fluxes. 

Alternatively, one could reduce the problem to a single variable, by considering linear 
combinations of the form 

*{x)=x#g + (l-x)*ffi Q . (69) 

For standard neutrinos and < x < 1 this combination should be positive. By requiring 
that the physical constraints are satisfied for any x one can determine the required proba- 
bilities/bounds. Both methods give the same results, which, taking into account the errors 
on neutrino cross sections, can be summarized as follows: 

• (a) the probability P\ for both and ^^NQ ^° ^ e positive is less than about 
2%. Should we disregard arbitrarily one of the experiments, one still has P\ < 6%, 7% 
or 9% neglecting, respectively, the results of Chlorine, Gallium or KAMIOKANDE 
experiment. This indicates that standard neutrinos ($3^ > and $cNO > ^) are 
unlikely; 

• (bl) to the 99.5% C.L., , without any a priori knowledge, + $cNO s fi° u ^ n °t 
exceed 0.7-10 9 cm~ 2 s _1 ; 
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• (b2) to the same confidence level, if one assumes a priori standard neutrinos, the 
combined flux of Be and CNO neutrinos does not exceed 2-10 9 cm -2 s _1 . 

Similar statements hold for 7 Be +CNO neutrino fluxes separately, see Table [TX| . 

The main message can roughly be summarized by saying that the probability of stan- 
dard neutrinos are low, not much more than 2%. However, the precise magni- 
tude of the probability should be taken with caution, at least for the following rea- 
sons: i) statistical and systematic errors have been combined together and ii) the assumption 
of Gaussian fluctuations might underestimate the probability, especially considering that we 
are dealing with the tail of the distribution. 



H. Experimental results and standard solar models 

Let us insist on the hypothesis of standard neutrinos and compare experimental infor- 
mation with theoretical estimates. 

We report in Fig. |7| the results of several recent solar model calculations [p!5H20|j23|j22]j56 



together with the experimental results. Some of the models predict a B neutrino flux close 
to the KAMIOKANDE value; however no model is capable of reproducing the low Be+CNO 
flux implied by the experiments. 

In Table we have considered only standard solar models [^,|T^,^,Q . For stan- 



dard neutrinos, the experimental information is also presented in the same Table. The 
discrepancy between theory and experiment is about a factor two for the boron flux. The 
discrepancy on $Bc+CNCH where the predicted values exceed the experimental 
upper bounds (99.5% C. L.) by a factor three, appears more important to us. 

The problem is mostly with beryllium neutrinos and let us examine it in some detail. 
The extraction of $g e from experimental data (with the requirement $CN0 — 0) yields an 
unphysically negative Be flux. Without any prior knowledge, $Be cannot exceed 1/10 of the 
RSM prediction at the 99.5% C.L. If we a priori force it to be non-negative, then the upper 
bound is 1/5 of the RSM at the 95% C.L.; a value as high as 1/3 of the RSM prediction 
is only allowed at the 99.5% C.L. All this indicates that 7 Be neutrino suppression is much 
stronger than that of 8 B neutrinos. 



I. The beryllium-boron relationship 

Additional insight on neutrino fluxes can be obtained by considering the physical con- 
nections among them. The relationship between between 7 Be and 8 B neutrinos, particularly 
emphasized by Berezinsky ||62|| , is most interesting. 

Both neutrinos are "daughters" of the 7 Be nuclei, see Fig. [8[ For this nuclide, electron 
capture (rate A e7 ) is favoured over proton capture (rate Ai 7 ), due to the absence of the 
Coulomb barrier. Thus the value of $g e is a clear indicator of the parent 7 Be concentration, 
n 7 : 

n 7 oc $ B eAe7 • ( 70 ) 
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Since $B e comes out to be reduced by a (large) factor with respect to the SSM prediction, 
the same reduction has to occur for the 7 Be equilibrium abundance (A e 7 is weakly dependent 
on temperature, and it is essentially scaled from measurements in the laboratory [|3*5|). The 
puzzle is thus with 8 B neutrinos, 



since: 
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(71) 



The observed (KAMIOKANDE) value of $3 being just a factor two below the SSM predic- 
tion, it looks that experiments are observing too high $g! To put it in another way, one 
cannot kill the father/mother before the baby is conceived. 

Should we insist on this approach, then we need to enhance Xn/X e j. Any attempt to 
reduce the discrepancy between the KAMIOKANDE and Chlorine experiments 
with respect to SSM by lowering Sn (the zero energy astrophysical factor for the 
p+ 7 Be — > 8 B +7 reaction) goes into the wrong direction. 

To make this argument more quantitative, let us define the reduction factors Ri with 



respect to the prediction of the Reference Solar Model {Ri = $j/$f SM ) 
Chlorine and KAMIOKANDE data one gets: 



R 



Be 



R 



< -1 ±0.8 



B 



From the 



(72) 



which is another way of presenting the 
KAMIOKANDE data, see Sec. [TTC 



'inconsistency" between the Chlorine and 



J. What if the sun were burning less now? 

One might speculate that the present luminosity L does not correspond to the present 
nuclear energy production rate L nuc in the sun. Actually it takes about eight minutes for 
neutrinos produced in the solar core to reach earth, whereas the time for electromagnetic 
energy to reach the solar surface is more than 10 4 years, and one might imagine that L nuc 
is different from L Q . Short time scale fluctuations in the nuclear energy production might 
not alter the photospheric temperature, the relevant time scale being of the order of 10 7 
yr because of the enormous amount of gravitational energy stored in the solar structure. 
Although this hypothesis of a thermal instability is rather extreme, let us consider it as a 
way of exploiting the full potential of the "solar model independent" approach. 

First of all, we remind that the result of combining Chlorine and KAMIOKANDE ex- 
periments (see Sect. |II Q ): 



~>Cl,int 



-0.48 ± 0.49 SNU (73) 



was independent of the solar luminosity constraint. The probability of finding Sci,int about 
lex below its physical limit is 16%. By exploiting Gallium result, we can find that the present 
situation is even less probable. If we insist that Gallium result is consistent with the other 
experiments, then the signal originates essentially from p (=pp+pep) neutrinos: 

$ p = = (60.2 ± 6.5)-10 9 cm" 2 s" 1 (74) 
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i.e. p neutrinos are essentially as many as calculated in standard solar models. From solar 
models we only assume $ pe p/$p = ^pep 1 /$p 5M , i.e. Q: 

$ pep = 2.36 ■ 10- 3 $ p . (75) 

One can thus calculate the contribution of pep neutrinos to the Chlorine signal: 

Sci, P e P = acW^pep = 0.23 ± 0.025 SNU. (76) 

After subtracting this contribution from Eq. ( [T3| ) we are left with: 

Sc^e+CAro = -0.71 ±0.49 SNU. (77) 

By requiring Sci,Be+CNO > 0, one immediately derives that the probability of the present 
situation is at most 8%. 

Essentially the same argument can be presented in a more precise form. We abandon 
the luminosity contraint Eq. (|3) or ( p4a| ) and still we have the three experimental signals, 



Eqs. (H - |6|) or Eqs. ( |54b| - |54d|) , in face of four unknowns: $ p , $ B e, $CNO an d 



One can thus express one flux (e.g. $Be) m terms of another flux (e.g. $cno) an d of the 
experimental signals. In this way one finds: 

$ Be + 1-77$CN0 = 4.59S CT - 1.29 • 10~ 2 S Ga - 5.04$ Bji , a 

= -3.0 ±2.2 , (78) 



which is essentially equivalent to Eq. (|77|) . This equation also shows the relevance of 
direct 7 Be neutrino detection. As an example, a measurement of $g e in excess of 3.3 
(10 9 cm" 2 s -1 ) would imply $CNO < at 3cr level, and thus would be a proof of non-standard 
neutrinos, even without the stationary sun hypothesis. 

K. Universal neutrino oscillations 

In the case of neutrino oscillations, similar arguments can be used when the averaged 
survival probability (P„ e _>„J is the same for all the components (pp, pep, Be, . . . ) of the 
neutrino flux. This situation is realized for Am 2 > 10 _3 eV 2 (in this case coherent matter 
effect are negligible at any point in the sun, for any component of the solar neutrino flux). It 
is important to observe that in this situation the cross sections <Ji averaged over the neutrino 
energy spectra are the same as for standard neutrinos. 

Thus by interpreting again $Be> ^CNO anc ^ as the electron neutrino fluxes at earth, 



the case of sterile neutrinos is exactly equivalent to the one just discussed in Sec. |II J| , 
and we get again Eq. (|7S|). This means that the probability of universal oscillations into 
sterile neutrinos is less than 8%. 

We remark that we used really minimal information from solar models, essentially Eq. 
([T5|). Clearly any additional hypothesis on the sun (e.g. a minimal non vanishing 7 Be flux) 



will essentially exclude this scenario, see also Ref. |63|] and Ref. [64|] 



ii 



We verified that £ varies by 2% when L is changed by 10%. 
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L. Concluding remarks 



All in all, a solar model independent evidence for non-standard neutrinos exists. It 
is however not overwhelming: the probability of standard neutrinos is less than a few 



percent (see Sec. II G), i.e. the indication is at the 2a level (we point out that this 



conclusion is reached when giving up all our understanding of stellar physics). 

• If one insists on standard neutrinos, then the 7 Be neutrino flux has to be drastically 
suppressed with respect to the prediction of SSMs and this suppression is stronger 
than that for 8 B neutrinos. Non-standard solar models attempting to solve the solar 
neutrino puzzle have to account for both these reductions. 

• All this does not imply that the Be flux on earth has to be small. That conclusion 
holds for standard neutrinos only. There are neutrino oscillation schemes which can 
account for all available data and, at the same time, predict a 7 Be signal quite con- 
sistent with the SSM prediction. (These models exploit the possibility of deforming 
energy spectra and/or transforming u e into u^, which are active in the KAMIOKANDE 
detector, see Ref. flB5H). In conclusion the direct detection of 7 Be neutrinos is 
crucial. 

In this respect, one has to avoid the temptations illustrated by the following story ||66 



The owner of a villa in Rome was explaining the civilization level of ancient Romans: "They 
even had telegraphs or telephones. When digging in my garden, I discovered ancient cop- 
per wires". His friend immediately went digging in his own garden and returned with the 
comment "You were right about their incredible civilization. They even had wireless com- 
munication. Indeed I found no cables upon digging" . 
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III. NON-STANDARD SOLAR MODELS: WHY AND HOW? 

A. Introduction 

For standard neutrinos, one finds that the fluxes of intermediate energy neutrinos 
( 7 Be and CNO) are strongly reduced with respect to the SSM expectations, so that the 
nuclear energy production chain appears strongly shifted towards the pp-I termination. 

The question addressed in the next sections is the following: is it possible to build non- 
standard solar models in agreement with available experimental data? In other words, if 
we insist on standard neutrinos, is the solar neutrino puzzle restricted to the results of SSMs, 
or is the problem more general? 

In order to enhance the pp-I termination, it is necessary that the ratio between the rates 
for the 3 He + 3 He and 3 He + 4 He reactions, 

R = WA^ (79) 

is drastically decreased with respect to the SSM prediction (here and in the following Ay is 
the fusion rate of two nuclides with mass numbers i and j and the average is meant over the 
energy production region). In the (realistic) approximation that the total flux is given by the 
sum of 7 Be and pp components, i. e. only ppl and ppll contribute, the first termination has 
probability P=A33/ (A33 + A34) and of course the other one has 1— P. In ppl two pp neutrinos 
are emitted, against one pp plus one 7 Be neutrino in ppll, so that: 

£bc __R_ 

$ PP 2 + R 

In the RSM one has = 0.09 so that Rrsm — 0.2; if the 7 Be neutrino flux has to be 

reduced by — say — a factor three, then one needs R ~ 0.07. 
One finds just two ways for decreasing R: 

a) Adjusting the parameters so as to lower the inner temperature. 

b) Adjusting the 3 He nuclear cross sections, so as to make the 3 He + 3 He reaction even more 
favoured with respect to the 3 He + 4 He reaction. 

We note, by the way, that the p+ 7 Be — » 8 B +7 cross section is not relevant for enhancing 
the pp-I branch. 

In this section we begin the discussion of case a), by identifying the parameters which 
could affect the inner solar temperature. Non-standard solar models with reduced central 
temperature will be explicitly presented in the next section, whereas the effect of varying 
the 3 He nuclear cross sections will be discussed in section [VJ As a way of characterizing the 
non-standard solar models which we shall discuss later, we introduce in the final part of this 
section some algorithms, essentially based on "homology" (scaling) concepts. 

B. Cooler solar models: why and how? 

By lowering the temperature and therefore the collision energies, the tunnelling proba- 
bilities are decreased for both the 3 He + 3 He and 3 He + 4 He branches, the latter being more 
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(80) 



suppressed as heavier nuclei are involved, see Eq. (|3|); in conclusion, the 3 He + 3 He branch 
gets favoured. 

How can one decrease the inner solar temperature? As we shall see, there are a few 
analogies between the solar core and the human body. For this latter the following statements 
clearly hold: 

i) the (absolute) temperature is fixed to the level of ± 1%; 

ii) alteration of the temperature is a symptom, and not an illness in itself; 

iii) once you measure the temperature somewhere, you know it everywhere. 

As regards the inner solar temperature, a comparison of the results of various SSM 
calculations (see Table |) immediately shows the first point. The stability of the internal 
temperature for any given input physics will be further analysed in Sec. [IV H| . Point iii) 



will be discussed more extensively at the end of this section. Let us now concentrate on 
ii). It means that one cannot treat — in principle — the central temperature T c as an 
independent parameter. In fact, T c cannot be decoupled from the solar structure, since 
severe constraints arise from the stability criterion. One has thus to study how T c is altered 
when physical/chemical parameters are varied, while the basic equilibrium conditions are still 
satisfied. In other words, one has to study "different-but-still-reasonable" suns, i. e. pseudo- 
solar structures where stellar matter behaves differently from the predictions of standard 
physics, the basic equations for stellar structure and evolution still being satisfied. 

In order to reduce T c one can resort to several manipulations, playing on the physical 
and chemical inputs which determine the structure of the star, see Sec. [TB| : 

a) A larger pp cross section. 
Increasing S pp (the astrophysical zero energy S-factor for the p + p — > d + e + + v e reaction) 
implies a lower central temperature, since fusion gets easier while the solar luminosity has 
to be kept constant. Although S pp is theoretically well determined (see Sec. |1V rj| ), one can 



introduce an artificial variation just to get cooler pseudo-suns. In the language of Sec. [IB 
one is essentially changing the energy production rate per unit mass. 
b) A less opaque sun. 

This is another way to get a cooler solar interior, as a smaller temperature gradient enlarges 
the region of nuclear burning, with less energy needed from the innermost core. In practice, 
this can be accomplished by using different ad-hoc assumptions: 

bl) the metal fraction Z/X is significantly smaller than that indicated by the photospheric 
and/or meteoritic composition; 

b2) the radiative opacity is smaller than that computed by several authors [BTipRBUJ ; 



bS) some new mechanism could contribute to energy transport through the sun. In this 
connection some years years ago (see e.g. p5| , pP|) the possibility of Weakly Interacting 
Massive Particles (WIMPS) captured by and trapped inside the sun was discussed. We 
do not consider this hypothesis in view of the negative results of direct searches for such 
particles and for the lack of any observational evidence in later stellar evolution |70|. 

According to Sec. [TB|, one could also study modification of the equation of state. We 
do not consider this possibility as in the region where neutrinos are produced the deviation 
from a perfect gas law, although relevant for a detailed evaluation of the solar structure, 
cannot deeply modify the present sun. 

On the other hand, there is an additional, still hypothetical possibility: 

c) A younger sun. 
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Should the sun be younger, the hydrogen mass fraction in the center would be higher and 
the same nuclear energy output could be produced at lower temperatures. 



C. Homology relationships. 



In the next sections we shall present many solar models, where some physical input X 
(solar age, chemical composition...) is varied from the starting value used in the calculations, 
X*, by a scale factor x: 



As a result, one finds profiles for the physical quantities Q {e.g. pressure, temperature, 
density...) which are different from the profiles Q* of the starting model, at any point in the 
solar interior. 

As we shall see, we find that in many cases the distributions of physical quantities follow 
with a good accuracy an "homologous" scaling relation 



(where m represents the mass coordinate m = M(r)/M , and M(r) is the mass within the 
current values of the radius r) i.e. for any physical quantity, the profile along the mass 
coordinate has the same shape as in the starting model; the difference with respect to it is 
just in a scale factor, which depends on the amount x of the variation of the parameter X 
and which can be different for different physical quantities (however fn = 1 for any Q and 
X, when X = X*). 

This is not in principle an unexpected behaviour, since we know from the theory of 
stellar structures that similar relations holds for "homologous" stellar model. To understand 
the point, let us recall that a stellar model is governed by a set of differential equations 
whose solutions give the correct distribution of the physical quantities Q throughout the 
structure. According to the linear structure of differential relations, given one solution 
(e.g. the starting model) further solutions can be generated by a scaling transformation of 
the physical quantities, as in Eq. (|82|) , the different scaling factors being related through 
algebraic constraints. 

Equation (|8~^) is thus essentially a way of transforming solutions into solutions. Note 
however that not any solution is obtained just by these transformations, much in the same 
way that not all stars are homologous to each other. As an example, we shall see that 
homology is violated when considering huge solar age variations, as the age alters the profile 
of helium abundance. 

More generally, bearing in mind that R Q has to be matched by tuning the mixing length, 
one finds that non-standard structures tend to have an homologous internal structure, with 
a strong departure from homology just in the more external layers. 

We recall that our main interest is in the resulting neutrino flux, so that two points have 
to be borne in mind: 

i) We are primarily concerned with the energy production region (m < 0.3 or R/R < 0.2). 

ii) We are mainly interested in the temperature profile, as neutrino production depends 
crucially on temperature. 



x = X/X* 



(81) 




(82) 
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All this means that we do not bother if Eqs. (|82|) are badly violated outside the central 
core, and that for physical quantities other than T we shall be satisfied if Eq. ( |82"D holds only 
to a fair approximation. What really does matter is the temperature profile in the energy 
production region. 

Clearly, if Eq. ( p2|) holds for the temperature profile, 

T(x,m) = f T (x)T*(m) , (83) 

then from the knowledge of the temperature at a point say the center, we are able to compute 
it at any other point, the only relevant parameter being the scaling factor: 

r = T c /T* c (84) 

In other words, the test of homology for the temperature profiles correspond to check the 
third statement mentioned in the previous section. 

It is useful to specify, in preparation for the next sections, the algorithms we use to test 
Eqs. (|82| ) and to extract the dependence of the physical quantities on the inputs which will 
be varied when building non-standard solar models. 

If Eq. (|82p holds for a physical quantity Q, then the ratio 

u ~ "TFtwT (85) 

is independent on the mass coordinate m and is purely determined by the physical input 
which has been varied. Qualitatively, this can be seen by looking at a graph where u is 
plotted versus m. For quantitative statements, it is useful to compute the average value of 
uj over the cells of our solar model and its variance Alv 

uj = -^—Y. u ( m k) ( 86 ) 



(^Y = —Y,Hm k )-ZJ] 2 (87) 

m k 

where we assume to divide (a portion of) the solar profile into N m cells, the fc-th one being 
centered at m = m^. By definition, for a perfect homology Auj = 0. The ratio 

5 n = Au/U (88) 

is an indicator of the validity of the homology relationship for the quantity Q, and typically 
we consider two regions: 

a) the energy production region (m < 0.3 or R/R < 0.2) which is of primary interest to 
us, as already remarked. 

b) The full radiative interior, (m < 0.98 or R/R Q < 0.7). Although such an extended area is 
not important for neutrino production, nevertheless it is interesting to study the behaviour 
of physical quantities up to the bottom of the convective layer, where useful constraints arise 
from helioseismological measurements. 
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D. Scaling laws for the physical quantities 



Once homology has been tested, one still has to study the function fn{x), i.e. the de- 
pendence on the input parameter X which is being varied. 

A natural parameterization, again reminiscent of those encountered in the study of the 
homology relationship is of the form: 

n(x,m) = x a n*(m) . (89) 

The coefficient a depends on the physical quantity Q (as well as on X). For simplicity of 
notation we will understand this dependence. Note that for x — > 1, Q — > Q*. 
From Eq. (|89]) , the power law coefficient a is determined as: 

0= '°f^> (90) 

log 2 

If we have built N x models, labelled by an index j specifying the value of the input 
parameter (x = Xj) and if each model contains N m cells, labelled by an index k indicating 
the value of the mass coordinate (m = m^) then one has N x xN m independent determinations 
of the coefficient a: 



log 



~ n*(m fc ) (gi) 

log Xj 

In order to extract a suitable average value for a, we proceed in the following way: 
a) for the j-th model, we perform an average over the cells 



1 



a 



(92) 



3 N 

m k 

and evaluate the corresponding variance A<x,- from: 

(A«,) 2 = ^E(^-S7) 2 - (93) 

m k 

b) We take then a weighted average over the different models, using l/(Aaj) 2 as a weighting 
factor: 

a ~ E,l/(Aa,) 2 (J j 

We remind that homology is well verified when Aaj/aj is small, and this provides a justifi- 
cation for the choice of the weighting factor. 

c) We can also define a variance Aa, specifying in some sense the global uncertainty on the 
coefficient just determined: 

(A«) 2 = ^E(« 3 -«) 2 (95) 
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Note that Aa/a essentially estimates the accuracy of the power law approximation (whereas 
Aaj/ctj indicates the accuracy of the homology relationships). 

Generally we consider two kinds of variations of the input parameter, for the calculation 
of the power law coefficients a and of their variances: 

i) small variations: | x — 1 |< 10%. This procedure, which was pioneered in |13[], is useful to 
study the effect of changing an input parameter of the SSM within its estimated uncertainty 
(which is generally of the order of few percent, see section P"V| ). 

ii) Large variations: typically by an order of magnitude, i.e. for x well outside the range 
allowed for the SSMs. It is in this way that one is really building and testing non-standard 
solar models. Such models are actually called for if one wants to effectively suppress 7 Be neu- 
trinos. 

Clearly the coefficients found in i) and ii) should be equal — - within numerical uncer- 
tainties — if the power laws were exact. Actually there is no deep reason for these laws 
to hold. They are just parameterizations of data and one should not be astonished to get 
different numbers for cases i) and ii). It is more of a surprise if the values are found to 
be close, indicating that the same simple parameterization holds over a wide range for the 
input parameter. 



E. Scaling laws for neutrino fluxes 

We assume again a power-law behaviour: 

$ . = x ai ■ $* 



(96) 



In order to determine the coefficients c^, we use algorithms similar to those of the previous 
section. As the fluxes are already summed over the cells we can skip point a) of the previous 
section and we construct directly the average over the models. Again omitting for simplicity 
the index specifying the flux component, one has: 



a 



AT 



N 



hg(xj 



(97) 



and the corresponding variance: 

(A«) 2 



N, 



E 



log 



log(x j 



a 



(98) 



Since all pseudo-suns have the same luminosity, they will give (approximately) the same 
total neutrino flux: 



J2Qi®i(x) = const = Y,Qi$* 



(99) 



so that equations of the form of Eq. cannot hold exactly for all the components of the 
neutrino flux and for arbitrary variations of the parameters. For small variations one gets, 
by differentiating Eq. (p^) with respect to x, the constraint: 
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i 



(100) 



which can be used as a check of the calculations. For the case of large variations it is 
convenient to use Eq. fl99|) as a way of expressing one of the fluxes in terms of the others, so 
as to maintain the luminosity constraint. This is best done for the case of pp+pep neutrinos, 
as these are the least sensitive to variations of the physical inputs. For these latter, thus, 
instead of Eq. fl^S|) we will generally use the expression: 

*P = % ~ T^Bc - $Bc) - %^($CNO - $CN0) (101) 



where $g e and $CNO are given through Eq. (|9 



F. Dependence on the central temperature 

All in all, we are mainly interested in the change of neutrino fluxes (and other physical 
quantities characterizing the stellar interior) on the inner solar temperature, when some 
input parameter is varied. 

This question can be easily answered on the basis of the above discussion, if we find - 
and we shall — that the temperature profiles satisfy the homology relationship: 

T(x,m) = '^-T*(m) . (102) 

c 

Neutrino fluxes, as well as the other quantities, can then be parameterized in terms of the 
scale factor T c /T*, see Eq. (jB4]): 

= $* . (r c /r;) ft (103) 
fi« = n* ■ (T c /T:f> . 

The Pi coefficients can be determined directly from the calculated fluxes for several 
models with a procedure similar to that presented in the previous section (this is the way 
we shall use for neutrino fluxes). Alternatively, one can profit from the previous results. If 
one has determined the dependence of the central temperature on the parameter x: 

T c = x aT ■ T* (104) 

and the dependence of the fluxes on the same parameter, see Eq. (|96l), then one has: 

rv ■ 

Pi = — (105) 

ax 

For neutrino fluxes, the most interesting comparison will be that of the values obtained 
by varying different physical inputs. In principle they do not need to be the same. If we find 
that they are close, no matter which parameter is varied, then this will be a confirmation 
of our expectation that T c is the quantity controlling the neutrino fluxes, independently of 
how that particular value of T c is achieved. 
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IV. LOW CENTRAL TEMPERATURE SOLAR MODELS 



A. Introduction 



Solar models builders claim that the central temperature of the sun is known with an 
accuracy of one percent or better. This claim is often questioned by other physicists, who feel 
such an accuracy as too high on a matter where no direct observational data are available. 
Independently of personal feelings, a few pertinent questions are the following. 

• If we insist that a low temperature solar model yields a drastically (say a factor three) 
reduced 7 Be neutrino flux, how much should the physical and/or chemical inputs of 
solar models be varied? Is it enough to go slightly beyond estimated uncertainties, or 
are wild changes actually needed? 

• Is it possible to get fluxes of both 7 Be and 8 B neutrinos consistent with available 
experimental information? 

Following the lines sketched in section [TII|, we now construct and discuss low inner tem- 
perature solar models. As the temperature is not an independent variable, we shall construct 
our pseudo-suns by acting on different inputs of the solar model: the p+p — > d + e + + v 
cross section, the metal content of the sun, the adopted values for the radiative opacity and 
the solar age. 

These inputs will be varied well beyond their estimated uncertainty, so as to build non- 
standard solar models. Generally we will attempt to vary the input parameters as long as 
we get a Be neutrino flux reduced to one third of the RSM prediction. 

To give and estimate of the sensitivity of central temperature to the chosen input param- 
eter, we will define x(0.1) as the value of the scaling factor such that the central temperature 
is reduced by 10%. 

A few common features of all the computed models will be summarized in section |1V Fj , 
namely: 

i) the temperature profiles appear, to quite a good approximation, homologous among the 
different models, in the sense specified in the foregoing section: 

T(m) = ^-T*(m) (106) 

ii) no matter how the temperature variation is obtained, the neutrino fluxes are essentially 
determined by the scale factor 

r = T c /T: (107) 

In the same section the temperature dependence of the main components of the neutrino 
fluxes, as obtained by numerical simulations, will be also demonstrated analytically. 

In section [IV G| we will speculate on the possibility of getting information on the central 
solar temperature through homology and of testing the homology relationship itself with next 



generation experiments, elaborating an observation by Bahcall |7T| , [72| , see also Ref . J73| . The 
final section contains our answer to the questions raised at the beginning. 
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In this entire section |V[ we shall for clarity divide the neutrino flux into the four com- 
ponents already introduced, $ p , $Be> ^CNO an d 3>b- I n ^ ne appendix B we will briefly 
discuss the ratios $pep/$p and ^n/^CNO mainly for substantiating our assertion about 
the stability of these quantities, among standard and non-standard models. 

Clearly our interest is on the changes of physical quantities with respect to the Reference 
Solar Model, BP95, when some input parameters are varied. 

Actually, following Refs. [71,[73] all throughout this section and section |V| we report 
results obtained by modifying the inputs of a starting solar model (CDF94) described in 
Ref. J73fl : the equation of state was taken from Ref. J75[| , internal opacity tables from Ref. [39|, 
corresponding to the chemical composition of Ref. [7B| and diffusion was neglected^. 

All quantities corresponding to the CDF94 model will be labelled with the 
index (*) here and in Sect. M. 

B. The p+p — > d + e + + v reaction rate 

The rate of the initial reaction in the pp chain is too low to be directly measured in 
the laboratory (even in the sun's center this rate is extremely slow, of the order of 10~ 10 
yr _1 consistently with the solar age) and it can be determined only by using the theory 
of low energy weak interactions, together with the measured properties of both the proton 
proton scattering and the deuteron. In terms of the astrophysical factor, Su (E) what really 
matters is its zero energy value, which for brevity, and following the usual notation, will be 
indicated simply as S pp . While we refer to Refs. 0,7^] for an updated review, we remark 



that the calculated values J77|,|7§ are all in the range (3.89-4.21) 10~ 25 MeV b, i.e. they 
differ from their mean by no more than 3%. Kamionkowski and Bahcall |77| give an useful 



parameterization, in terms of the three quantities of physical interest for the determination 
of S pp : the squared overlap integral A 2 (0), the ratio G^/Gy of the axial to vector coupling 
constants and the fractional correction 5 to the nuclear matrix element due to exchange 
currents 

S pp [l(T 25 MeV b] = 3.89[A 2 (0)/6.92][(GVGV)/1.2573] 2 [(1 + 5)/lM} 2 . (108) 

The most recent evaluation of A 2 (0) is from J77[] , obtained by using improved data for pp 
scattering and for the deuteron wave function and also including the effect of vacuum polar- 
ization. The estimated uncertainty is about ±1%. The ratio Ga/Gv can be obtained with 
an accuracy of about 0.3% from a weighted average of five precise modern experiments. 
The contribution of the exchange currents is 5 ~1%, with a comparable uncertainty. 



In summary, one has J77| 



S pp = 3.89 • 10^ 5 (1 ± 0.01) MeV b . (109) 

Although some warning is in order as to the meaning of the quoted (la) error, one 
may conclude that well known physics determines S pp to the level of few per cent or even 



12 CDF94 cannot be considered anymore as a standard solar model, since it does not satisfy the 
helioseismological constraints Eqs. (||) and (|8|). 
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better. Nevertheless, as explained in the previous section, the variation of S pp well beyond 
its estimated uncertainty provides a good theoretical laboratory for investigating alternative 
solar-like structures. 

In the RSM as well as in CDF94, Spp corresponds to the central value of Eq. (|109|) . 

When varying S P p(Spp — ► s P pSpp), we considered both the case of small variations (s P p 
in the range 0.9 to 1.1 ) and large variations, up to s pp = 3.5 which corresponds to a 
7 Be neutrino flux reduced by a factor three. For drastically reductions of $b, unreasonable 
variations of S pp are needed, orders of magnitude larger than compatible with the estimated 
uncertainty. Briefly, we remark the following occurrences (see Ref. J7D| for details): 

i) in the energy production region, temperature, density, pressure and radius all satisfy the 
homology relationship Eq. (|88|) to better than 1%, and the same holds throughout all the 
radiative interior (see Fig. |9|). The same holds over the entire explored range in s pp . 

ii) The hydrogen mass fraction, as a function af the mass coordinate, is essentially unchanged 
with respect to the CDF94 estimate, as it is constrained by the solar age. 

iii) For the power law coefficients a and j3 [see Eqs. (^) and ( |103|) 1 of the quantities charac- 
terizing the physical interior we found the values in Table |X|. The power laws look accurate 
(Aa/a and A/3//3 being just a few percent) and the coefficients for small and large variation 
of Spp are consistent, within their estimated uncertainties. In addition s pp (0.1)=2.5, i.e. one 
has to multiply S pp by this huge factor to reduce the central temperature by 10%. 



iv) The dependence of neutrino fluxes is also shown in Table [X[ For the case of small 



variations, already investigated by Bahcall in Ref. [|l^,^J,^y], we essentially agree with his 
results. Even for large variations, power laws are very accurate for $g e and $g. For $ pp the 
accuracy is smaller, for the reasons outlined in section [T.I fB[ This holds for $CN0 as well, 
as this flux is the sum of two terms ($n and $q) which depend differently on temperature. 

We stress that these features (points iii and iv ) can be well understood analytically, 
assuming that the starting model and the pseudo-suns are connected by a homology trans- 
formation and that the hydrogen mass fraction profiles are the same, see Ref. 



C. Radiative opacity 



An extensive and critical discussion of the uncertainties on radiative opacities is given 
in Ref. |14|]. In the energy production region the typical difference between the outputs of 
the Livermore and the Los Alamos code is about (2-5)% |8T|. At least half of the opacity in 
the central region is due to scattering on electrons and inverse bremsstrahlung in the field 
of H and He nuclei, processes which can be calculated with an accuracy of about 10%, or 

rjj estimate a la uncertainty of about 2.5%. On the 
claim the uncertainty to exceed 5%, a point criticized 
argued that some plasma physics effect have not 



better. Bahcall and Pinsonneault 
other hand Turck-Chieze et al. 182 



m 



14f . Recently Tsytovitch et al. |83 



been included in the calculations of the Livermore group, so that opacity at the solar center 
might be overestimated by 9%. 

Although it is hard to make a definitive statement on such a complex matter, we shall 
conservatively take 5% as a lcr uncertainty. 

Since we aim to lower the internal temperatures, we investigated the effect of reducing 
the opacity. We scaled it uniformly along the solar profile by a factor opa with respect to 
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CDF94. We consider opa as low as 0.6. This corresponds to a temperature reduction of 
6% and $b g ~ 0.5<I>g e . For even smaller opacity, the resulting pseudo-sun would have an 
original helium abundance well below the cosmological value. By extrapolating, one finds 
ppo(0.1)=0.42. 

The following points are to be noted: 
i) The homology relationship is accurate to at least 1% for any variable characterizing the 
internal structure, over the entire radiative region, see Fig [TT]. ii) Density, pressure and 
radius are essentially insensitive to opacity variations, see Table [XI]. On the other hand, 
when the opacity decreases the hydrogen mass fraction increases. This can be understood 
by observing that as the star gets less opaque, the interior becomes cooler. At the lower 
temperature, the pressure gradient needed to sustain gravity is then maintained with a larger 
hydrogen abundance. 

iii) The dependence of neutrino fluxes on the opacity parameter and the connection with 
the central temperature are also shown in Table [Xl. Again our results for small variations 



are well in agreement with those of [13,24,80 



D. The "metallicity" of the sun 

Let us remind that the metallicity^ is an input parameter constrained by photospheric 
observations. For a comprehensive review about the metal content of the solar interior we 



refer again to |15] and we mention just the main points: spectroscopic observations give the 
mass abundance (Z) of the heavy elements, relative to the hydrogen mass fraction (X), in the 
atmosphere of the present sun. About 75% of the heavy elements is accounted by carbon, 
nitrogen and oxygen. For these elements, the analysis of vibration-rotation and pure rotation 
lines of molecules as CO, CH, OH, NH in the infrared, from space experiments P^ , [76]J55|] 
provide accurate information. The abundance of elements heavier than oxygen is usually 
determined by absorption lines in the optical range; these clearly show the predominance of 
Iron. The relative distribution of these elements can be usefully tested to the distribution 
of CI carbonaceous chondrite meteorites, which should keep the composition of these not 



volatile elements in the original solar nebula [HRfEJ. It is important to note that the recent 



photospheric Iron abundances |45| agree now well with the meteoritic values. 



According to the most recent evaluations j|5] , one has for the photosphere of the present 
sun: 

(Z/X) phato = 0.0245 (110) 

with an accuracy better than 10%. 

We remark that, due to diffusion towards the solar center, the original heavy elements 
abundance in the sun should be higher, by 10-15%. 

When building solar models, this time we keep z = (Z/X)/ (Z/X)* as a free variable. We 
considered small variations (z = 0.9-1.1) and large variations, with z as small as 0.1. With 



13 According to the astrophysical jargon, any element heavier than He is termed "metal" or "heavy 
element" 
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such a small value the central temperature is decreased by about 10% and $g e is about 1/3 



of the initial prediction. In this case z(0.1)=0.17. While we refer for details to [73,74]], we 
summarize here the main results, see Fig. [0]. 

i) The temperature profile satisfies the homology relationship with an accuracy better than 
1% throughout the entire radiative interior. Homology holds also for R(m), but with less 
accuracy. On the other hand, even in the energy production region, homology is only a fair 
approximation for pressure and density: for large variation of Z/X even in this restricted 
region 5uj/oj 5%. 

ii) The picture is somehow similar to that of the opacity variations: R remains essentially 
unchanged on the average, whereas X grows as metallicity is diminished, see Table [XI 1| . 

iii) This similarity also reflects on the neutrino fluxes, the temperature dependence of these 
latter are essentially the same as in the previous subsection a part from $CN0> see below. 

As already mentioned, the similarity with the variations of the opacity can be readily 
understood since metallicity affects mainly the stellar opacity. 



E. Rejuvenated suns 

Unlike the other parameters used to constrain the SSM, the age of the sun is not an 
observable. It is inferred from the dating of the oldest meteorites, provided that a connection 
between the formation time of the meteorites and the birth of the sun (i.e. the ignition of H 
burning) is achieved. A recent discussion is provided by Wasserburg in ||15|| . On this basis 
Bahcall and Pinsonneault estimate (at la) 

t = 4.57 ±0.01 Gyr (111) 

Since the main difficulty is in establishing the evolution phase of the sun at the formation 
time of meteorites, we shall take as a conservative estimate of the uncertainty the duration 
of the pre-main sequence phase, At Q ~ 30 Myr. 

The solar luminosity and age fix the amount of H which is burned into He in the solar 
core. In younger suns more H is available, so that fusion reactions are more likely and the 
observed luminosity can be reached at smaller T c . 

Correspondingly, a younger sun has a different composition in the energy production 
region, this difference becoming more and more marked as t is shortened, and this drives 
the structure progressively away from homology (see Fig. O). 

In CDF94 we started with t* = 4.6 Gyr. 

We reduced the solar age down to t Q = 0.1 t , this extreme corresponding to a central 
temperature reduction of about 6% and $b c « 1/2 $g e . For even shorter ages, the structure 
of the pseudo-sun would be deeply modified by the occurrence of a central ( 3 He driven) 
convective core. We found the following results: 

i) homology is now just a fair approximation (see Fig. [13|). Even restricting oneself only 
to the energy production region, the accuracy in temperature is about 2% or less; that in 
density and pressure about 10% and for the hydrogen abundance it is merely about 25%. 

ii) Correspondingly, the power laws for the variables characterizing the interior are only 
rather approximate (A/3//3 ~ 20%, see Table Kill]) . 
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iii) As regards neutrino fluxes, we are in fair agreement with the results given in |T3|j2^j80[] for 
small variations. The dependences on temperature are in any case similar to those previously 
found. 



F. Low temperature models and neutrino fluxes: summary and explanation 

An important general feature of the models discussed above is the approximate homology 
of the temperature profiles 

T{m) = rT*(m) , (112) 

where m = M/M is the mass coordinate, and the factor r depends on the parameter which 
is varied but does not on m. 

We have verified that Eq. ( |1 12j ) holds to an accuracy better than 1% in the entire radiative 
zone (M/Mq < 0.98 or R/R < 0.7) for all the models considered, except for huge (and 
really unreasonable) variations of t (see Fig. [14]). It is worth noting that T(m)/T*(m) stays 
constant throughout a region where T(m) change by a factor five. The scaling factor r may 
be taken as a the ratio of the central temperature T c to that of the starting model: 

t = TJT* c . (113) 

The coefficients /3j for the power laws of neutrino fluxes vs temperature 



are collected in Table |XV| , where we include for completeness all flux components. 

One notes that (3 p , /3g e , and (3b are largely independent of the parameter which is being 



varied (see also Fig. [15]). In other words, these fluxes are mainly determined by the 
central temperature, almost independently of the way the temperature variation 
is imposed. 

The dependences of the fluxes can be understood semi-quantitatively by simple analytical 
arguments. As a zeroth order approximation, let us assume that energy production occurs 
entirely through the ppl termination, and that the chain is fully equilibrated. Requiring 
that the rate of 3 He burning (nl(av) 33 ) corresponds to the fusion rate {n\{av)n/2) at each 
point in the stellar core, the equilibrium 3 He density n 3 is given by |79fl : 

2 {crv) 3 3 Qpp{av) 33 

where e is the energy production rate per unit mass. 

The production of 7 Be nuclei, and consequently of 7 Be neutrinos, can be treated as a 
perturbation to the ppl termination. Practically every 7 Be nucleus produced is destroyed 
through electron capture, with emission of a 7 Be neutrino. The production rate per unit 
mass of the latter, W7 Be , is thus equal to the production rate of 7 Be nuclei, again per unit 
mass: 
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W7 Be = n 3 n 4 (av) u /p (116) 
With Eq. ( |115| ) and n 4 = pYN& one gets: 

W'Be = T= , =^V^P ( 117 ) 

4 \/<5pp V(0-V>33 



The nuclear reaction rates ((tu)^ are strongly temperature dependent. This dependence is 
usually parameterized by power laws [BiJ : 

(av) i:j ocT^ (118) 



and the coefficients 7^ are given in Table |XIV . 

Clearly quantities like Y, e and p are also connected with temperature, but the depen- 
dence is much weaker than for the reaction rates; as a further approximation, we assume 
they are the same as in the starting model. On the other hand, we use homology for the 
temperature profiles. In this way, from Eq. (|1 1 7| ) , the neutrino production at any mass 
coordinate m is related to that of the starting model: 



1 , 



wr Be (m) = w? Be (m)r 734 -2 733 (119) 
The same equation obviously holds for the fluxes on earth, $j = — ^ / dmwi(m): 

ES 

$Be = $Be^ 4 " 1733 = *Be^ (120) 

One can study the production of 8 B neutrinos similarly. Their production rate per unit 
mass is: 

ws B = n 1 n 7 (av)i 7 /p , (121) 
whereas for 7 Be neutrinos one has: 

wr Be = n e n 7 (av) e7 /p . (122) 
By eliminating the equilibrium 7 Be nuclei density n 7 one has: 

(<™)i7 ni , . 

W»B = ~, v W7 Be (123) 

{(TV) e7 n e 

Again assuming n\jn e to be essentially that of the starting model, one obtains after inte- 
grating over the mass coordinate: 

We are now also able to estimate the temperature dependence of $ p . We recall that 
the two main components are $ p and $Be; an d that their sum is fixed by the luminosity 
constraint. By differentiating with respect to T c , this implies approximatively: 

Pp = -/^% = -°- 6 ( 125 ) 

The values we find are in agreement with the numerical estimates of Table |XV| . It is not 
surprising that the analytical values are quite close to the coefficients obtained by varying 
S pp , since scaling works best in this case. 

Let us remark a few relevant points: 
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the temperature dependences of $ p ,$B e and $g are well under control. One sees 
from the foregoing discussion that they are essentially determined by the behaviour 
of (crv)ij as a function of temperature; the latter behaviour is fixed mainly by the 
Coulomb barrier ||35|| . Solar physics only enters through the homology relationships of 



the temperature profiles. 

The flux $B e can been determined independently of the value of the 7 Be lifetime. No 
matter what the value of this latter is, practically all 7 Be nuclei produced will emit a 
7 Be neutrino. 

The ratio of $g to $B e is essentially governed by nuclear physics, its temperature 
dependence being that of (av )n/(crv ) e7 . 



For any temperature one has 



r (126) 



i.e. the suppression/enhancement of $g is much stronger than that of $Bc- 
A reduction of $B e to 1/3 of the RSM prediction implies a reduction of 
$B by an order of magnitude. This essentially illustrates the failure of low 
temperature models when compared with experimental data. 



The CNO flux appears sensitive not only to the temperature, but also to other parameters 
characterizing the solar interior. The drastically different exponents, found when varying 
the metallicity, can be understood by noting that the efficiency of the CN cycle is affected 
not only by temperature but also by the number densities of nuclei which act as catalyzers 
of the chain. 



G. More about homology: additional consequences and a possible test 



One generally thinks that the neutrino production zone is so well hidden below the solar 
surface that it can hardly be studied experimentally other than with neutrinos. 

Actually, the accurate homology of temperature profiles, valid up to the border of the 
radiative interior, indicates a strict connection between the properties of the energy produc- 
tion region and of more external layers of the sun. If one is confident in homology, then a 
measurement of temperature at, say, the bottom of the convective zone immediately gives 
the temperature of the solar center. 

A new generation of experiments is being planned for detecting monochromatic neutri- 
nos produced in electron capture ( 7 Be +e~ — > 7 Li+z/) and in the pep (p + + p — > d + v) 
reactions 0,^,^]. Bahcall [|7T] , [72|| has pointed out that one can, from measurements of 
the average energy difference between neutrinos emitted in solar and laboratory decay, infer 
the temperature of the production zone. The possibility of measuring inner solar temper- 
atures through thermal effects on monochromatic neutrino lines is extremely fascinating 
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(although remote). In this respect the homology relationship, Eq. (|112|) , is particularly in- 
teresting, see Fig. |16[ If homology holds, then a measurement of the solar temperature in 
the 7 Be production zone gives the value of T c . 

In addition, the homology relation itself is testable - in principle - by comparing the 
temperatures at two different points, as can be done by looking at the shapes of both the 
vse and Upep lines. We remark that this would be a test of the energy transport mechanism 
in the inner sun. 



H. On the accuracy of the central solar temperature 

A rough estimate of T c can be obtained by equating the thermal energy of a hydrogen 
nucleus to its gravitational energy {kT c « GM Q m p /R Q ). In this way one finds T c « 2 • 10 7 
K, in good agreement with the much more refined SSM estimates (see Table |). Different 
standard solar models give the same value of T c within 1% and solar model builders claim 
that the present accuracy on T c is of the same order. 

From the preceding discussion, the reader can derive his own opinion about this claim. 
We have seen that the main regulators of the central solar temperature are S pp , the metal- 
licity, the solar age and the opacity of the solar interior. In terms of the uncertainties in 
these quantities, one has: 

AT C 2 ASpp 2 AOpa 2 A(Z/X) 2 AAge 2 

(-jr) = (°*A^) +(0*,^— ) (127) 

Using the power law coefficients olt% and the estimated uncertainties summarized in Table 



XVI|, one actually gets AT C /T C »s 1%. One also sees that most of the error results from the 



uncertainties in opacity and in metallicity (note that with the uncertainties as estimated in 



T4| , [T5|| the global uncertainty on T c would be halved). 

The next question is how the main parameters should be varied in order to get a solar 
model with drastically reduced 8 B and/or 7 Be neutrino fluxes. For instance if one requires 
$b/^b ~ 1/2 a reduction of T c by 3% is needed. Although this does not look terribly 
outside the allowed T c range, really huge variations of the physical inputs are needed (see 
again Table |XV1| , 4th column) . 

The situation is even more desperate if one tries to reduce the 7 Be flux by a factor 
two (three). In this case temperature has to be reduced by at least 7% (13%), which 
requires the really unreasonable variations given in the last columns of Table |XVI] . The 
input parameters are to be varied by an order of magnitude (or more) with respect to their 
estimated uncertainties. Even a 40% reduction of opacity can only reduce the 7 Be neutrino 
flux by at most a factor two. Last but non least, as is clear from Eq. (|126|) , it is essentially 
impossible to account for the reductions of both $g e and $g. 

One concludes that a solution to the solar neutrino problem cannot be found 
merely by reducing the central temperature . 
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V. NUCLEAR REACTIONS IN THE SUN 



One may suspect that the solar neutrino problem is due to some inadequacy in our 
understanding of nuclear reactions in the solar interior. What could be wrong with the 
nuclear burning rates used in standard solar model calculations? There are (at least) three 
sides to this question, viz.: 

• Nuclear physics: as repeatedly stated, the astrophysical S-factors used in stellar model 
calculations are generally obtained by extrapolating experimental data taken at ener- 
gies higher than those relevant for the solar interior. Although the underlying theory 
is robust, one can be suspicious of extrapolations. 

• Atomic /molecular physics: experiments in the laboratory use atomic or molecular 
targets. At the lowest measurable energies, electron screening is relevant and its effect 
has to be subtracted when deriving cross sections for bare nuclei. Indeed, the effect 



of elctron screening has been detected [S8|, however theory and experiments seem to 
disagree |88|j89|] . 



Plasma physics: the burning rates for bare nuclei are then to be corrected for the 
screening of nuclear charges by the solar plasma. In any calculation, the predicted 
effects are small; however, the theory is not completely satisfactory [9(J and there are 



no direct experimental tests. The disagreement between theory and experiment for 
electron screening in atomic/molecular targets provides some warnings: although that 
is a different context, one has to keep an open mind about the possible effect of plasma 
screening on neutrino production. 

As shown in Sec. f|, the main problem is with the 7 Be neutrinos, so that particular atten- 
tion has to be given to the reactions in the pp chain preceding the formation of 7 Be nuclei; in 
other words, those reactions which determine the branching between the ppl and ppll chains 
(see Fig. [I]). The role of the p+p cross section was already discussed in Sec. [IV B| , as it is one 
of central temperature regulators. The value of the p+d — > 3 He +7 cross section is unimpor- 
tant. It is orders of magnitude larger than that of the weak interaction p+p^ d+e + + v e 
process, so that on a very short time scale equilibrium is reached and any practically pp 
reaction is followed by the production of one 3 He nucleus. 

The cross sections for reactions between He isotopes, 

3 He + 3 He -> 4 He + 2p and 3 He + 4 He -> 7 Be + 7 , (128) 

are clearly crucial and we shall pay particular attention to them. On the other hand, the 
proton capture by 7 Be : 

p + 7 Be ^ 8 B +7 (129) 

is of minor relevance for the flux of 7 Be neutrinos and governs essentially only the production 
of 8 B neutrinos, see Sec. |V D\ Similarly, the raction: 

p+ 14 N^ 15 + 7 , (130) 
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which will be briefly discussed in Sect. [VE| , essentially determines the production of CNO 
neutrinos only. 

The goal of this section is to study the changes induced in neutrino fluxes when the 
nuclear burning rates are changed with respect to RSM inputs. The comparison of these 
results with neutrino experiments will be discussed in the section [VT]. As mentioned in 



Sec. IV, we recall that our starting standard solar model here is CDF94 174,73 



A. The status of 3 He + 3 He — > 4 He + 2p reaction 



The relevant energy range is determined by the energy E Q and the halfwidth, A/2 of the 
Gamow peak^. For the central solar temperature one has E Q =22 keV and A=12 keV. The 



available experimental data pl|-p7i1 are shown in Fig. [T7|. Data below 25 keV do not exist 
and at the lowest measured energies errors are of the order of 20% or greater, so that some 
extrapolation is necessary in order to reach the relevant energy range. 

The astrophysical S-factors are usually parameterized with a Taylor expansion: 



s l3 (E) = s^o) + s' tl (o)E + -s>;mE> 



131) 



where mainly the coefficient 6^(0) matters for the sun, as E Q is much smaller than the 
nuclear energy scale. In Table [XV1J] we present some estimates of 6 33 (0): the value used 
in the Caughlan and Fowler compilation jy] corresponds to the experimental result of [92 
Parker and Rolfs |9B| give a weighted average of several experimental results. 



We derived a new estimate of 633(0) by reanalyzing all available data, except for those 
of the pioneer experiment by Good et al. |JTJ which is systematically a factor 2-3 below 
the others and carries no estimated error. A fit to the data using Eq. ( p. 31]) gives 6(0) = 
(5.3±0.1) MeV b. The data at low energies are however affected by electron screening ||35|| . 
This effect is negligible at E > 100 keV and data above this threshold give the value of the 



4th row in Table [XVII| . A slightly smaller result is obtained if one considers all energies 
and corrects for electron screening effect in the adiabatic approximation as shown in 
the last row. We consider this last value as the best estimate of 633(0) for bare nuclei. 
Correspondingly we find: 



S 33 (E) = (5.1 ± 0.2) + (-3.0 ± 0A)E + ^(3.0 ± 1.0)E 2 . 
This expression (energies in MeV and 6 in MeV b) gives a good fit to all data (Xd.o.f. 



(132) 



0.8) 



At the Gamow peak near the solar center, one thus finds 633 
than, but still consistent with the value used in the RSM, 633 
used 63*3 = 4.98 MeV b. 



5.0 MeV b, slightly higher 
4.8 MeV b. In CDF94 we 



14 



E = \.22{Z\ ZtfiTg) 1 / 3 keV and A = 0.749(Z 1 2 Z^T^) 1 / 6 keV, see e.g. || 
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B. The status of 3 He + 4 He — > 7 Be + 7 reaction 



The range of relevant energies is essentially the same as above (E Q = 23 keV, A = 12 
keV), whereas data are available only at E > 100 keV. Concerning the experimental results 
collected in Fig. 



lq , the following comments are needed: i) the original data of |99[ have 
been corrected, following the comment in ||100|| ; ii) data from ||101|| have been multiplied 
by a factor 1.4 according to [102]; iii) errors quoted in [103] have been doubled so that 



2|; iii) errors quoted in 

fluctuations among data points become statistically consistent; iv)at all the energies where 
data are available electron screening is irrelevant. 



Again what matters is the astrophysical S-factor at zero energy. In Table |XV111| we report 



some different determinations of this quantity. The value used in ||T| coincides with that 
of the review paper 
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which was obtained as a weighted average of the extrapolations 
provided by different experiments. A similar, more recent analysis []98| gave a slightly smaller 
value. It has however to be remarked that the extrapolations were performed using different 
theoretical models, so that combination of extrapolated values is dubious. We performed 
a new analysis of all experimental data ||1 05| , p9Kl 1| , |10 3| Jl 6|J1 2| , |i~07| , |l 8| | . For a quadratic 
expansion of the astrophysical S-factor we find 



S 34 (E) = [(4.8 ± 0.1) + (-2.9 ± 0.2)E + (0.9 ± 0.1)E 2 ] ■ 10" 



(133) 



(again E in MeV and S in MeV b). Alternatively, by using an exponential parameterization, 
as frequently adopted in the literature and supported by theoretical models (see e.g. 
we get: 



S U {E) = (5.1 ±0.1) ■ 10- 4 exp[(-0.83±0.07)£ + (0.25 ± 0.03)£ 2 



(134) 



By using these two parametrizations, at the Gamow peak near the solar center one finds 
respectively S 34 (E ) = 4.74 and 5.01 in 10" 4 MeV b. One sees that uncertainties due to the 
extrapolation procedure are at least comparable to the quoted statistical error, so that the 
global error is about ±2 ■ 10 -5 MeV b. 

The value of Su(E ) used in the RSM is based on J9~8"| , after introducing a 1.6% decrease 
due to the vacuum polarization effect ||110|1 , yielding Ss 4 (E ) = 5.17 • 10~ 4 MeV b, slightly 
larger than that given by our preferred expression Eq. (|134|) , but still consistent with it 
within uncertainties. In CDF94 we had S^(E ) = 5.26 • 10" 4 MeV b. 



C. Neutrino fluxes and helium reactions 

In order to estimate the dependence of 7 Be and 8 B neutrino fluxes on the nuclear cross 
sections, essentially we repeat here the argument of Sec. [IV R Let us first consider the local 
equilibrium concentration of the parent 7 Be nuclei, n 7 . 7 Be is created in the 3 He + 4 He reac- 
tion and destroyed, essentially via electron capture, with a lifetime r e7 . Thus at equilibrium: 

n 7 = r e7 n 3 n 4: (av) 34 , (135) 

where n« is the number density of the nuclei with mass number equal to i. The 3 He equi- 
librium density is obtained by equating its creation rate (i.e. the rate of the p+d-^ 3 He +7 
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reaction, which equals the rate of p+p— > d +e + + z/ e ) to the burning rate, which is dominated 
by the 3 He + 3 He reaction. With this approximation, one gets: 

By using the above equations one can derive: 

n 7 = —7=T e7 nin A \l {av) n i (137) 

If the rates (crv)^ and (0^)34 differ from the starting inputs, one expects that only 
the probabilities of the pp terminations are varied whereas the densities n 1; n 4 and the 
temperature are essentially unchanged with respect to the starting predictions, since the 
hydrogen burning rate and the helium abundance in the solar interior are determined by 
the present luminosity and the present age of the sun. This expectation is confirmed by 



numerical experiments, see for example Fig. [19]. One thus gets that only the following 
combination of (crv) 33 and (av)^ matters: 

X = 4^ (138) 



f (<™)33 

In terms of x one has: 

n 7 = n** (139) 
X 

We remark that Eq. ( |139|) holds at each point in the solar interior. 

Assuming for the moment that there are no resonances in the energy range of interest 
and that the low-energy astrophysical S-factors differ from those used in the starting model 
by a constant factor s^- (Sij(E) = SijS*j(E)), then at any point in the sun: 

X 534 (140) 



X V S 33 

and consequently: 



nj = n ssMj^ (141) 

V^33 



The 7 Be and 8 B neutrino fluxes, which are obviously both proportional to n 7 , scale then in 
the same way 



$Be = $Be s 34 8$* ; $ B = $B *34 s£* (142) 



The p (= pp+pep) neutrino flux can be best derived by the luminosity constraint, Eq. (|13|). 

Numerical experiments confirm these analytical estimates. In Fig. [H] one sees that 
for each component the flux is actually determined by the variable x an d that Eq. (|142|) 



provides a good approximation. In Table |XiX| we show the power law coefficients for the 
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main components of neutrino flux. Our results for small variations agree with those which 



can be derived from ||24j| . The last column presents the case of large variations. One notes 
that the numerical values are well in agreement with the analytical estimates presented 
above. 

Note also that pp+pep and CNO neutrinos are essentially insensitive to variations of \ 
as expected since the temperature is unchanged. 

Large variations are really necessary if one wants a reduction of 8 B and 7 Be neutrino flux 
by a factor 2-3. As an example, if S33 is kept fixed then £34 has to be reduced by a factor 
2-3, or S33 has to be enhanced by a factor 4-9 if S34 is unchanged. This is clearly in conflict 
with the experimental situation discussed in the previous sections. 

The only way out, a desperate one, is to invoke a resonance [|1 1 1|] in the 3 He + 3 He reaction, 
the resonance energy E r being below the experimentally explored region, i.e. E r < 25 keV. 
Such a resonance is not predicted theoretically; furthermore experimental searches for excited 
6 Be states in reactions like 6 Li(p,n) 6 Be and others failed []3"5]| . One can however not definitely 
exclude this possibility [j35|l . 

A resonance would affect the various components of the neutrino flux differently. Qual- 
itatively, a very low energy resonance will be more effective in the more external (cooler) 
solar regions, so that the 7 Be neutrino flux can be more suppressed than the 8 B flux. The 
converse is true for a higher energy resonance, the border between the two regimes being 
provided by the Gamow energy for 3 He + 3 He reaction, about 22 keV near the solar center. 



The best case for our purpose [|112|| is that of E r =0 which corresponds, for a strong 
to 



$Be/$ 



Be 



0.75 (143) 



B 



In practice, the effect of such a zero energy resonance can be mimicked by introducing a 
resonance S-factor, S res and parameterizing the fluxes as: 



$ B = (l + |?) (145) 

All other fluxes are unchanged, except for the pp+pep flux which can be derived by the 
luminosity constraint, Eq. Q4"3|). In the Sec. |V]]we compare the results of solar models with 
such a hypothetical resonance with experiments. 

Before concluding this section, we remark that, in principle, also a narrow resonance 
below threshold (i.e. E r < 0) could work, since again it is more efficient in the 7 Be production 
region than in the solar center, where 8 B is produced. Note however that the resonance has 
to be very close to the threshold, otherwise it is uninfluential. 

The experimental situation on such hypothetical resonances should be clarified by the 
ongoing LUNA experiment |113| at the underground Gran Sasso National Laboratory. 
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D. The p+ 7 Be -> 8 B +7 reaction 



The determination of Sn, the astrophysical factor for the p+ 7 Be — > 8 B +7 reaction, 
involves several complications: i) extrapolations are needed to reach the relevant energies in 
the solar interior (E = 19keV) as measurements have been performed only at Eqm > 100 
keV ||114[ - p.l9|| , see Johnson et al. ||120|| for a detailed discussion about the extrapolation 



procedure; ii) the number of 7 Be nuclei present in the target is most accurately determined 
by monitoring the build-up of 7 Li (the 7 Be decay product) as a function of time, using the 
7 Li(d,p) 8 Li reaction, see ||121]1 - For this technique to be useful, the cross section for this 
latter reaction must be known and there has been considerable work on this point since 
1978, resulting in a determination with an accuracy of 6% [ 121|| . Note that in the past 
different values of this normalization cross section have been used by different authors. 

The available data from different experiments, all normalized to the same value a dp = 
157 ±10 mb and extrapolated according to ||120| , are summarized in Table [XX| . For deriving 
an average value, the experimental results of ||116|| and ||119|| are particularly important since 
they correspond to measurements performed over a wide energy interval, reaching the lowest 
energies. One has to remark that, at each energy, the values of ||119|| are systematically lower 



than those of [[L16|| . Keeping into account this uncertainty, Johnson et al. derived^: 



S 17 (0) = 22.4 ± 2.1 eV b 



(146) 



This is the value adopted in the Reference Solar Model (Table |XXV111| ). 

The Coulomb dissociation process has recently attracted a great deal of attention as 
an alternative method to study radiative capture reactions of astrophysical interest at low 
energies. The process can be treated as the absorption of a virtual photon, essentially 
the inverse of the radioactive capture process, see Fig.[H|. The cross section for Coulomb 
dissociation of 8 B — the 208 Pb( 8 B , 7 Be p) 208 P reaction — was measured with a radioactive 
8 B beam at RIKEN [|124j| and a preliminary value 6*17(0) = 16.7±3.2 eV b was deduced from 
the data. A theoretical analysis of the same data by Langanke and Shoppa | ]125[ | yielded an 
even smaller value, Si7(0) = 12 ± 3 eV b, which was however criticized in [|126| . 

All in all, this indirect approach looks very nice, but the extraction of Sn is experimen- 
tally difficult and theoretically complex, and we agree with the authors of Ref. 1125] that 



11 The recently developed technique of Coulomb dissociation might prove itself as a very useful 
tool in nuclear astrophysics . . . The 8 B Coulomb dissociation experiment at RIKEN takes a 
first step in this direction and should be continued and refined. However a reliable deter- 
mination of the astrophysically important cross section for the 7 Be (p^f'B reaction from 
8 B Coulomb dissociation experiments has to wait until improved data become available" ||127| . 

We thus consider Eq. ( |146| ) as our reference value, and we will study the effects of varying 
S17. It is clear that acting on a very minor termination of the fusion chain, the properties of 
the solar interior (profiles of temperature, density, . . . ) will be unchanged. This expectation 
is confirmed by numerical experiments, see for instance Fig. 22. Thus, only the relative 



15 Recently, the value o^p = 146 mb has been recommended in Ref. [ |122| . A reanalysis of all data 
by U. Greife and M. Junker within the NACRE collaboration yields S 17 (0) = 19.9 eV b pi. 
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intensity of the ppll and ppIII terminations are affected, and consequently only the 7 Be and 
8 B neutrino fluxes will change, their sum being fixed at the SSM value. If the astrophysical 
factor is scaled by an amount sn, one has: 

^ " -ir£ (147) 



*B + $Be = $ B + $ Be 



This clearly implies: 



$ R = <5>S SM Sl7 1+# B /# Be ( U8 ) 
^B ^B 6l 7 I+S17** /** V 1 * ) 
' B' Be 

1+** /<6* 

$ R = &SSM + B Be 
^Be ^Be l+Si7*g/* Be 

We remind that the ppIII termination is very disfavored with respect to the ppll, so that, 
to a very good approximation: 

$ B = s 17 $| ( 149 ) 
$ Be = $ Be 

These analytical estimates are (obviously) confirmed by numerical experiments, see Ta- 
ble [X-X-lj . Note that the slight dependences of $ pp and $CNO are n °t significant, being at 
the level of the numerical accuracy. 

In conclusion, only the flux of 8 B neutrinos is significantly affected when 5x7 
is changed, the dependence being linear to a very good approximation. Clearly, 
playing with Sn cannot be the solution of the solar neutrino puzzle, as the 
production of intermediate energy neutrinos is unchanged. 



E. The p+ 14 N -^ ib O +7 reaction 



We will briefly review the status of this reaction since it is the slowest process in the 
main CN cycle, and thus it controls the production of CNO neutrinos. 

The Gamow peak in the solar center corresponds to E D = 27 keV, whereas laboratory 
measurements have been performed at energy in excess of 100 keV, so that extrapolations 
are necessary. 

The data of earlier investigations |128| - |133|1 were extrapolated by various groups, yield- 
ing different values for S , i i 4 (0) ranging from about 2 to 10 keV b This discrepancy 
illustrates the difficulty in determining absolute cross sections. Fowler et al. in 1975 || 1 3 1|| 
recommended S^i^O) = 3.32 keV b, by means of some unspecified, judicious treatment of 



experimental data and/or theoretical input. Schroder et al. ||135|| performed a comprehensive 
experimental study, the only one covering continuously the energy range E = 200 — 3600 keV 
and measuring absolute cross sections, 7-ray angular distribution and excitation functions. 
This work removed most of the appearent discrepancies among previous experiments. The 
extrapolated result, Si^O) = 3.20 ±0.54 keV b, is " essentially identical with the previously 
recommended value of 3.32 keV b" |p8| , |135|| . In his latest compilation Fowler adopted 



the value found by Schroder et al. On the other hand, Bahcall and Pinsonneault |TJ| quoted: 
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51,14(0) = 3.32 ± 0.40 keV b (150) 
and later reduced the central value by 1%, after including the effect of vacuum polarization 

n • 

Concerning solar neutrinos, the CNO neutrino flux scales linearly with 51,14: 

$CNO = $ CNO s i,i4 ( 151 ) 
In order to keep the same luminosity, also pp and pep neutrino fluxes are slightly changed: 

*p = *l a vF > ( 152 ) 

all other components being essentially unsensitive to S1.14. 

F. Plasma screening of nuclear charges 

1. The results of different models 

The burning rates of bare nuclei derived from experiments need to be corrected to take 
into account the screening provided by the stellar plasma. The study of screened nuclear 
reaction rates was started with the pioneer work of Salpeter ||136|| ; it has been addressed 



by several authors, see e.g. [ 137 -140], and recently reviewed in [ |141|]90|1 . In the sun the 



screening effects are small, however different calculations yield relatively different nuclear 
reaction rates. 

As a starting point let us neglect any screening effect, i.e. reactions take place for bare 



ions with rates Xnos- The results of the corresponding solar model are shown in Table [XXII 
(NOS). Due to the screening, the actual rates A will be larger : 

A = AW (153) 

where the enhancement factors / depend on the reaction and on the plasma properties. 
Various approaches have been developed for evaluating these factors. 



In the weak-screening approximation (WES), originally introduced by Salpeter ||136|| , one 
has for a Debye plasma wherein partial electron degeneracy is included: 



\nf w * b = Z x Z 2 e l j{R D kT) (154) 

where Zi j2 are the charges of the reacting nuclei, T is the temperature and Rn is the Debye 
radius, see ref. ||142|| . In this scheme, the reacting particles are assumed to move slowly in 



comparison to the plasma particles (adiabatic approximation). Also modifications of the 
Coulomb potential are assumed to be sufficiently weak so that the linear approximation 
holds (see [|120|| for a more extensive discussion about the validity of this approach). In 



the sun, the weak-screening approximation is justified (to some extent) for the pp-reaction, 
whereas the other nuclear reactions occur in the so-called intermediate screening regime. 
Graboske et al. |I3§ (GDGC) used Eq. (g5§ when f WES < 1.1; for larger values (up 



to 2) they derived the enhancement factors by using general thermodynamic arguments and 



interpolation of Monte Carlo calculations. The explicit expressions can be found in ref. |138 . 
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Mitler |140|| (MIT) developed an analytical method which goes beyond the linearized 



approach and which correctly reproduces both the limits of weak and strong screening. 
Neglecting the small effects of a radial dependence in the effective potential, see [141|, the 
enhancement factors are given then by: 

ln/ M7T = -\^en e fR% [(d + C2 + 1) 5/3 - (Ci + 1) 5/3 - (Ca + l) 5/3 ]/(^T) (155) 

where Ci,2 = ^Zi^/ {^n e B? D ) and n e is the average electron density 

In the Salpeter approach one assumes the reacting nuclei to be so slow that the plasma can 
fully rearrange itself while the nuclei are moving (as in the Born-Oppenheimer approximation 
used in molecular physics) and dynamical effects (corrections to the Born-Oppenheimer 
approximation) are completely neglected. Carraro et al. |139| (CSK) observed that the 
reacting nuclei actually move faster than most of the plasma ions (the Gamow peak energy 
is generally larger than thermal energy), so that ionic screening plays a smaller role under 
this condition. They calculated the dynamic response of the plasma in the framework of the 
linearized theory. The resulting enhancement factors are expressed in terms of those of the 
weak screening: 

lnf CSK = C lnf WES . (156) 

The coefficients C essentially specify the corrections to the adiabatic limit, which obviously 
corresponds to C — 1. At the center of the sun the correction factors are ||139|| : C p+P = 0.76, 



C'iHe+ 3 He = 0.75, Cs He+ 4. He = 0.76, C p+ 7Be = 0.80, C p+ i4 N = 0.82. Note that the isotopic 
dependence is rather weak (in all previous models there was no isotopic dependence). 

We shall not discuss the strong screening limit which is definitely too far from solar 
conditions (see e.g. IP-42H ). 

From Table [XXI 1| , where we report the results of solar models corresponding to the 
different approaches, one notes the following features: 

i) The largest differences arise between the no screening model and the weak screening 
model. $g can vary by at most 15%, the Chlorine signal is stable to within 13% and the 
Gallium signal at the level of 3%. 

ii) The GDGC model, extensively used in stellar evolution codes, yields values very close 
to the no screening model; the difference between the two is at the level of 1% for $g e and 
$B, as well as for the Chlorine and Gallium signals. 

In Fig. ^ we show the enhancement factors along the solar profile, calculated by using 
the different prescriptions outlined above, for the reactions relevant to hydrogen burning in 
the sun; as concerns the CN cycle, we pay attention only to the slowest reaction: p + 14 N 
^ 15 0+7. All the enhancement factors depend very weakly on the mass coordinate, at least 
in the energy production region (M/M Q < 0.3). This is clear in the weak screening regime, 
since the dependence on the solar structure parameters is just of the form p/T 3 (see Eq. 



( |154j) and remind Rd oc JT/p) and this quantity is approximately constant along the solar 
profile. The same holds in the strong regime ||142j| , and thus the approximate constancy 
in the intermediate regime is not a surprise. For these reasons, in the following we shall 
concentrate on the enhancement factors calculated at the solar center (see Table [XXI 11| ). 



The weak-screening approximation, Eq. ( |154| ), always yields the largest enhancement 
factors, as it is physically clear since electrons and ions are assumed to be free and capable 
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of following the reacting nuclei and in addition the electron cloud is allowed to strongly 
condense around the nuclei (in the linear approximation the electron density becomes infinite 
at the nuclear site). By using the Mitler model, where electron density at the nuclear site is 
fixed at n e , one obtains smaller enhancement factors. The same holds for the model where 
the limited mobility of ions and thus their partial screening capability is taken into account. 
The GDGC enhancement factors are systematically smaller than the others (except for the 
pp reactions where, by definition, they are equal to the weak screening prescription). It is 
thus clear that the corresponding neutrino fluxes and experimental signals are the closest 
ones to those of the no-screening models. One notes that the enhancement factors for He+He 
and p+ 7 Be are very close: actually, in the weak screening approximation only the product 
of nuclear charges enters. 

All in all, the enhancement factors are relatively close to unity, however none of the 
approaches to screening discussed above is completely satisfactory. The weak screening 
approximation is not justified for reactions other than the pp, since ZiZ 2 e 2 / (RpkT) is 
not small. The GDGC results stems from an interpolation of numerical computations and 
the prescription of the authors yields an unphysical discontinuity at the border between 
the weak and intermediate regimes ||141|| . The CSK result, which incorporates dynamic 



effects of finite nuclear velocity, is however derived in the framework of a linear theory, 
i.e. the weak screening approximation. The Mitler approach goes beyond the weak screening 
approximation; on the other hand, the partial mobility of ions due to ion-interaction effects 
and/or to the finite thermal velocity is not taken into account. Also, the value of the electron 
density at the nucleus is somehow artificially kept equal to the average electron density n e . 



2. A model independent analysis 
The role of screening on solar neutrinos can be investigated, more generally, in a model 



independent way [p0| . In the previous section we saw that a few features are common to 
any screening model: 

i) the enhancement factors / can be taken constant in the energy production region, so 
that one needs to specify the values at the solar center only. 

ii) They are almost insensitive to isotopic effects, i.e. approximately 

/ 3 He+ 3 He = / 3 He+ 4 He ( 157 ) 

and we will refer generically to an enhancement factor for helium-helium reactions, /He+He- 

iii) To a good approximation, the enhancement factor is determined by the product of 
the electric charges of the reacting nuclei, so that we can take: 

/ p + 7 Bc = /He+He ( 158 ) 

In this case one is left with just three numbers, / p + p , /He+He an d / p + 14 N' which we 
shall consider as free parameters. We recall that / > 1 since in the plasma the Coulomb 
repulsion between the reacting nuclei is diminished. 

The introduction of (spatially constant) enhancement factors is equivalent to an overall 
change of the astrophysical S-factors: 
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Sij > Sij fi+j (159) 

and we can exploit the results for the variations of the astrophysical S-factors presented 
previously. 

Concerning the role of /p+p, we recall that an increase of S pp immediately implies a 
reduction of the central temperature T c . From Table [X|: 

T c /T c NOS = (/p+p)- 1 / 9 (160) 

where the superscript NOS refers to the no-screening solar model. From the same table, 
by using the dependence of neutrino fluxes on Spp, one immediately gets the effect of the 
screening factor: 

*Be = $Be° 5 (/ P + P r L1 (161) 
$B = K° S (/p+p)- 2 ' 7 (162) 

$CNO = $CNO (/ P + P r 2 - 2 ■ (163) 

Concerning the 3 He + 3 He and 3 He + 4 He reactions, we recall from the foregoing section 
that the 7 Be equilibrium concentration scales as: 

n 7 oc S u / y/s^ oc yj / He +He (164) 

This is clearly the dependence of 7 Be neutrino flux; for $g one has to remind that it is 
proportional to Syj and thus an extra power of /ne+He occurs: 

®Bc = Ke S (/Hc+Hc) 1/2 (165) 

*B = *g° 5 (/He+He) 3/2 (166) 
The enhancement factor for p+ 14 N only matters for the CN cycle: 

*CNO = *CNO / P + 14 N (167) 
We can put together all the previous results in the following way: 

*Be = *B? S (/He+He) 1/2 (/p+p)" 1 " 1 (168) 
®B = <° S (/He + He) 3/2 (W^ 7 (169) 

*CNO = $ CNO /p+ 14 N (/p+p)" 2 ' 2 ■ (I 70 ) 

The behaviour of pp-neutrinos, as usual, can be best derived best by using the conser- 
vation of luminosity. 

By using the above equations with the enhancement factors given in Table [XX111| , one 
can quantitatively reproduce, to a large extent, the numerical results presented in Table 
KXTj. Furthermore, we have verified that these analytical results are quite accurate for a 



wide range of the / factors, by using our stellar evolution code. 

Note also that although 3 He + 3 He and 3 He + 4 He have the same enhancement factor, 
the equilibrium concentration of 7 Be, and thus $g e and are changed when screening is 
introduced. 



63 



VI. NON-STANDARD SOLAR MODELS AND EXPERIMENTAL RESULTS 



A. Introduction 



We compare now the predictions of non-standard solar models with the experimental 
results on solar neutrinos. The basic questions are: 

• Is there a solar model that accounts for all available results, assuming standard 
neutrinos? 

• What would change if one of the experiments were wrong? 

• If — as it happens — no model is successful, what is the reason of the failure? 



For this purpose, we consider solar models characterized by three parameters: 

(1) the central solar temperature, T c , which accounts for (most of) the effect of changing 
the astrophysical factor S pp for the p+p— > 2 H+e + +i/ e cross section, the solar opacity or the 
age of the sun; neutrino fluxes are determined essentially by T c , independently of the way 
in which that particular temperature is achieved (see Fig. [T5| and Refs. ( |73| , |143| , |TP| ); 

(2) the astrophysical factor S33, which can be used as an effective parameter controlling 
both cross sections for the He+He reactions: by varying S 33 , at fixed S34 , one can tune the 
parameter x — which determines the neutrino fluxes when the cross sections for 
the 4 He + 3 He and/or 3 He + 3 He channels are altered (see Sec. [V C|) ; 

(3) the astrophysical factor Sn for the p+ 7 Be — > 8 B + 7 reaction. 

These three parameters are really independent of each other; we recall that T c is essen- 
tially unaffected by variations of S33 and/or S17. 

Possible variations of screening factors will not be discussed extensively since, as shown in 
the previous section, they can be rephrased in terms of variations of astrophysical S-factors 
(see however Ref. p0|). 



The dependence of the neutrino fluxes on the three parameters is shown in Fig. For 
T c we use the exponents corresponding to (large) variations of S pp (see Table [X]). For the 
dependence on the cross sections of the two He + He reactions, we explicitly consider the 
case of a zero-energy resonance in the 3 He + 3 He channel, as given by Eqs. ( |144| ) and ( |145| ) 
and use the resonant contribution S^ s as a free parameter. This possibility has the best 
chance of producing agreement with the data, since it suppresses more strongly $g e than 
$g: a nonresonant variation is even less effective ||111| , |112| , [73| . The pp- neutrino flux, <3> pp , is 



determined by imposing the luminosity constraint, see Eq. ( |101| ). Note that S17 affects only 
the 8 B neutrino flux significantly. As reference fluxes $f SM , we use the ones of the reference 
standard model (RSM) of Bahcall and Pinsonneault ("best model with diffusion"), which 
are given in Table [Hj (BP95). 

We remark that the precise values of the exponents and coefficients in Fig. ^4| are not 
important for our discussion; somewhat different exponents could origin, if the change in T c 
was induced by means other than changing the astrophysical factor S pp , or if these power 
laws were fitted on different ranges of parameters. Such different choices do not affect the 
essence of our conclusions. 
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In the following sections, we first vary one single parameter at a time, and then analyze 
all possible combined variations. In addition, we repeat the analysis after having arbitrarily 
discarded any one of the experiments and, finally, we discuss the reasons of the failure of all 



our attempts. This section updates and extends the analyses of Refs. 



73 


74 


145 


147 



B. Fitting all experimental results 

Let us consider the ensemble of solar models originating fluxes $j as parameterized in 



Fig. 24. We performed a \ 2 analysis to establish quantitatively how well these models 
compare with experimental data. A given solar model predicts theoretical signals {X = 
Gallium, Chlorine and KAMIOKANDE) according to the formula 

Sf = J2a Xii ^ , (171) 

i 

with the averaged neutrino cross sections <Jx,i given in Table 0. These theoretical signals 
are to be compared with the experimental signals reported in Table IVT. GALLEX and 
SAGE results have been combined into a single Gallium signal as in Eq. fl26|) . As usual, we 
define the likelihood function as: 

X 2 (T C , SS-, S„) = - S*)V£(S? - S?) • (172) 

XY 

The covariance matrix Vxy takes into account both experimental and theoretical uncertain- 
ties. These latter always include the errors on the averaged neutrino cross sections <Jx,i- 
Theoretical uncertainties on neutrino fluxes reflect uncertainties on T c , S^ s and Sn. When 
some of these variables are used as free parameters only the uncertainties corresponding to 
the remaining ones are included. The propagation of these latter uncertainties to the fluxes 
is established by means of power laws similar to those of Fig. optimized for the case of 
small variations. The matrix Vxy is not diagonal because the same parameter can affect 
more than one flux and the same flux contributes in general to more than one signal. The use 
of the full covariance matrix is really necessary, since otherwise apparently good fits can be 
achieved in an unphysical way |73|J1 44|JT48|J1 49|| . Error correlation means, for instance, that 



we cannot use the uncertainty in $g to strongly reduce its contribution to the Davis experi- 
ment, while having at the same time a smaller reduction in the KAMIOKANDE experiment. 
More details on the calculation of the covariance matrix can be found in Ref. [f73j|. 

We explored a wide region for the three parameters: T c down to 8% of the RSM value, 
Si7 from zero up to five times the RSM value and S^ s from zero up to 200 times the value 
of S% 3 SM . With three experimental results, one has two degrees of freedom (d.o.f.) when 
the three parameters are varied one at the time, and no d.o.f. when all three are varied 
simultaneously. The results of this analysis are summarized in Table [X-X1V| which needs the 
following comments. 

(1) None of the attempts succeeds in giving an acceptable x 2 even when we use as many 
parameters as the experimental data. 

(2) The smallest \ 2 ■, obviously obtained when all parameters are left free, corresponds 
anyhow to quite unphysical variations: T c comes out 7% smaller than the RSM value and 
the Siy astrophysical factor is 5 times larger. 
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(3) All the best fits give too high signals for the Gallium and Chlorine experiments and 
too low a 8 B flux as compared to KAMIOKANDE; in other words, the KAMIOKANDE 
signal looks too high for an astrophysical solution. 

(4) When varying only one parameter at a time, a resonant increase of the 3 He + 3 He cross 
section is the most effective way, as it allows a higher suppression of 7 Be relative to 8 B neu- 
trino flux. The central temperature is less effective, since it gives a too large suppression 
of the 8 B neutrino flux. S17 only affects the 8 B neutrino flux and it is totally useless for 
the present solar neutrino problem,which is mostly a problem of the intermediate 
energy neutrinos. 

(5) A significantly reduced \ 2 is obtained by changing S i7 together with T c or S 33 . Note, 
however, that the smallest \ 2 is obtained at the border of the parameter space, i.e. in an 
extremely unphysical region. On the other hand, there is no significant gain when varying 
at the same time T c and Sg| s as they both mainly shift the balance towards the ppl chain. 



C. What if one experiment were wrong? 

There is no rational reason for doubting any of the experimental data. Nevertheless in 
the same spirit as discussed in Sec. [TI], we have repeated the analysis of the previous section 
excluding in turn the KAMIOKANDE, Chlorine and Gallium result. (Note that, excluding 
the Gallium result, we are giving up two experiments, GALLEX and SAGE.) 



The results of this exercise, reported in Table |XXV| , show that the situation is essentially 
unchanged: 

(1) Again none of the attempts succeeds in giving an acceptable x 2 , even when we use 
more parameters than experimental data. 

(2) As previously, the best fit corresponds anyhow to extremely unphysical values of the 
input parameters. 

(3) The smallest x 2 are obtained disregarding either the KAMIOKANDE or Gallium 
data. In other words, these two results taken together are the hardest to reconcile with 
astrophysics. The point is that Gallium result implies that only the ppl chain 
is effective, whereas KAMIOKANDE shows that ppIII, and hence ppll, are 
operational. 



D. What went wrong with solar models? 



Tables XXTV and XX V| present the quantitative evidence that any solar model (which 
we are able to parametrize) cannot account for the experimental data, even when one of 
these data is arbitrarily removed. In this section we discuss the physical motivations of 



such failure, giving through Figs. ER [26] and 127] a graphical illustration of the results of 



Tables XXIV and XXV 



The starting point is that any combination of two experimental data requires, with 
respect to the RSM, a strong suppression of the intermediate-energy-neutrino flux ( 7 Be and 
CNO), much stronger than that of the 8 B flux, as was extensively discussed in Sec. [TT]. 

In the ($b> ^Be+CNo) plane, assuming standard neutrinos, the regions allowed at the 
2a level by each experiment and their intersection are indicated in Fig. |25. 
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The aim of non-standard solar models is to push the predictions close to this area allowed 
by experiments. 

In Fig. |25] we show the effect of changing one of the three parameters in turn. A reduction 
of Si 7 only affects the 8 B flux and we are left with a much too high $Be+CNO Reducing the 
temperature decreases both $g and $Be; however $g is suppressed more strongly than $g c 
contrary to experimental evidence. The resonant increase of the 3 He + 3 He cross section 
succeeds in reducing $b 6 almost to zero while leaving part of $g; however, $Be+CNO i s still 
too high since $CN0 = ^CNO * s itself sufficient to spoil the agreement with the data. 

In Fig. |2B| we show the effect of changing two parameters at the same time. The best 
case is when T c and S i7 are changed. In fact, only T c is able to cut both $b 6 and $CN0; 
and Si 7 helps to bring $g back up. We stopped at T c = 0.88 and S17 = 5Si 7 SM (which are 
already completely unphysical values!) and yet have not reached the allowed region (which 
of course could be attained by even smaller T c together with larger S17). 

Finally, we show the effect of changing all three parameters in Fig. ^j]. This case does 
not differ from the previous one; the main difference is that since we use S|f s to cut $g e and 
use T c mostly to cut •I'cNO we can reach even a lower value of $Bc+CNO with a relatively 
higher temperature. However, within the framework of present physical knowledge, one 
cannot find a justification to accept a non-standard solar model with the required decrease 
of temperature (7%). 

For the same reasons, playing with screening factors does not help; see Fig. ^ where the 
three screening factors / p + p , /He+He an d / p + 14 n ( see section [VF 2| ) are arbitrarily varied 
by huge amounts with respect to the RSM estimates. Again the suppression of intermediate 
energy neutrinos is too weak in comparison with experimental results. 

All in all, the fits we found are not acceptable, i.e. the chances of an astrophysical 
solution look very weak even if we are extremely generous on the region spanned by the 
physical inputs. In Sec. |I| we found that the probability of standard neutrinos is at most 
2% when all fluxes are left as absolutely free parameters. It is thus not a surprise to reach 
an even stronger conclusion when some astrophysical input is used. 
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VII. THE NEVER ENDING STORY? 



The sun or the neutrino, who is at fault? To express our opinion on such a matter, let 
us summarize the main points of this review. 

Assuming standard neutrinos: 

• the available experimental results appear mutually inconsistent, even if some of the 
four experiments was wrong. This is a solar-model independent evidence for non- 
standard neutrinos. This evidence, however, is not an overwhelming one, since there 
is still a few percent (at most) probability of neutrinos being standard, i.e. we have a 
"2-<T level effect" (see Sec. |I|). We remark that this conclusion is reached at the cost 
of giving up all our understanding of stellar physics. 

• Even neglecting these inconsistencies, the flux of intermediate energy neutrinos ( 7 Be + 
CNO), as derived from experiments, is significantly smaller than the predictions of 
SSMs. The main puzzle is with 7 Be neutrinos, for which the theoretical predictions 
are really robust. If one insists on standard neutrinos, one has definitely to abandon 
Standard Solar Models, all predicting a much too high 7 Be neutrino flux. 

• In addition, the suppression of 7 Be neutrinos is much stronger than that of 8 B neutri- 
nos. Non-standard solar models fail to account for both these reductions. 

• Even extremely non-standard models cannot account for data, even when one of the 
experimental data is arbitrarily omitted. All the fits that we attempted are unaccept- 
able (see Sec. |VT| ), i.e. the chances of a nuclear or astrophysical solution look to us very 
weak, although we have been extremely generous in the region spanned by the physical 
inputs. Since we have found that the probability of standard neutrinos is about 2% 
or less when all fluxes are left as free parameters, it is not a surprise to reach an even 
stronger limit when some astrophysical input is used. 

Let us wait (and hope) for future experiments, bringing direct, decisive evidence of some 
non-standard neutrino property. 
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APPENDIX A: OUR STANDARD SOLAR MODEL 



We present the solar model resulting from an updated version of FRANEC^ (Frascati 
Raphson Newton Evolutionary Code), where helium and heavy elements diffusion is included 
and the OPAL equation of state (EOS) |36|J150|| is used. The EOS is consistent with the 



adopted opacity tables, i.e. the most recent evaluation from the same Livermore group ||151 



and this should further enhance the reliability of the model. In addition, updated values 
of the relevant nuclear cross sections are used and more refined values of the solar constant 
and age are adopted. 

We discuss the effects of each of these improvements, showing that they are essential 
in order to get agreement with the helioseismological information about the bottom of the 
convective zone. We also calculate neutrino fluxes and the expected signals in ongoing 
experiments. 

FRANEC has been described in previous papers (e.g. see Ref. | 152 , 153|| ). Recent determi- 
nations of the the solar luminosity (L Q = 3.844- 10 33 erg/s) and of the solar age (t Q = 4.57T0 9 
yr) are used ||15|| . The present ratio of the solar metallicity to solar hydrogen abundance by 
mass corresponds to the most recent value of Ref. j45| : (Z/X) p ^ oto = 0.0245. 

Following the standard procedure, for each set of assumed physical inputs, the initial Y, Z 
and the mixing length parameter a were varied until the radius, luminosity and (Z/X) photo at 
the solar age matched the observed values within a tenth of percent or better. We considered 
the following steps: 

a) As a starting point we used the Straniero equation of state UTEfl , the OPAL opacity 
tables which were available in 1993 |4^|j39[1 for the solar metallicity ratio of Ref. [45], combined 
with the molecular opacities of Ref. [ |154j| ; diffusion was ignored. This model may be useful 
for a comparison with the Bahcall and Pinsonneault model without diffusion, described in 
ref. as the chemical composition is the same, although it is not the most updated one. 

b) Next, we introduced the OPAL equation of state. With respect to other commonly 
used EOS, this one avoids an ad hoc treatment of the pressure ionization and it provides 
a systematic expansion in the Coulomb coupling parameter that includes various quantum 
effects generally not included in other computations (see refs. |36j and ||150|| for more details). 

c) We used the latest OPAL opacity tables [ 151 1 , solar metallicity ratio as in pS| and 
again the molecular opacities of Ref. ||154j| for temperature below 10 4 K. With respect to 



Ref. |40| the new OPAL tables include the effects on the opacity of seven additional elements 
and some minor physics changes; moreover the temperature grid has been made denser. 

d) We included the diffusion of helium and heavy elements. The diffusion coefficients 
have been calculated using the subroutine developed by Thoul, see Ref. |^TJ. The diffusion 
equations were integrated numerically. At any time step, after updating the physical and 
chemical quantities as usually done in FRANEC, we took into account the effect of diffusion. 
Our time steps are about 3 • 10 7 yr. The variations of the abundance of H, He , C, N, O 
and Fe are followed; all these elements are treated as fully ionized. According to Ref. all 



16 



Based on work of Ciacio et al. 
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other elements are assumed to diffuse at the same rate as the fully-ionized iron. To account 
for the effect of heavy element diffusion on the opacity coefficients, we calculated the total 
heavy elements abundance in each spherical shell in which the solar model is divided and 
we interpolated (by a cubic spline interpolation) between opacity tables with different total 
metallicity (Z = 0.01, 0.02, 0.03, 0.04). 

e) Finally, we investigated the effect of updating the nuclear cross sections for 3 He+ 3 He 
and 3 He+ 4 He reactions, following a recent new analysis of all available data, described in 
Sec. |V} For S34, we used the exponential parametrization of Eq. (|134j ). 

The resulting solar models are summarized in Table |XXV1| , and deserve the following 
comments. 

(a — > b): The introduction of the new OPAL EOS reduces appreciably the initial helium 
abundance. The Straniero (1988) EOS underestimates the Coulomb effects neglecting the 
contribution due to the electrons, which are considered as completely degenerate, whereas 
the OPAL EOS includes corrections for Coulomb forces which are correctly treated. The 
models with Straniero EOS have a higher central pressure and a higher central temperature, 
and correspondingly a higher initial helium abundance. The effect of an underestimated 
Coulomb correction was discussed in Ref. ||17|| . Note that the transition between radiative 



and convective regions is not correctly predicted by the model, the convective region being 
definitely too shallow. 

(b — ► c): The updating of the radiative opacity coefficients has minor effects. The 
convective zone is again too shallow. 

(c — > d): This step shows the effects of diffusion. Helium and heavy elements sink relative 
to hydrogen in the radiative interior of the star because of the combined effect of gravitational 
settling and of thermal diffusion. This increases the molecular weight in the core and thus the 
central temperature raises. The surface abundances of hydrogen, helium and heavy elements 
are appreciably affected by diffusion. For example, the initial value Y in = 0.269 is reduced 
to the present photospheric value Y p hoto = 0.238, The predicted depth of the convective zone 
and the sound speed are now in good agreement with helioseismological values. It is natural 
that when diffusion is included the radiative region reduces its size. With respect to models 
without diffusion, in the external regions the present helium fraction is reduced while the 
metal fraction stays at the observed photospheric value. Thus the opacity increases and 
convection starts deeper in the sun. 

(d — > e): The modifications of our solar model arising from the new values of the nuclear 
cross sections are negligible with respect to the other improvements just presented. 

The predicted neutrino fluxes and signals are summarized in Table [XXVI I| for our models. 
All in all, the results are quite stable with respect to the changes we have introduced as 
long as diffusion is neglected. On the other hand, due to the higher central temperature, 
model d) has significantly higher 8 B and CNO neutrino fluxes. It is essentially the increase 
of $b which enhances the predicted Chlorine signal. The slight change in the nuclear cross 
sections weakly affect neutrino fluxes and signals. The 7 Be and 8 B fluxes are reduced by 
about 5% as a consequence of the correspondingly smaller value of S 34 (see Table |XXV1J| ). 
Should we use the polynomial expansion of Eq. ( |133| ) one would get a further 5% decrease. 
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APPENDIX B: THE STABILITY OF PEP AND CNO NEUTRINO FLUXES 



We already compared the predictions of Standard Solar Models by different authors in 
Sec. H see also Table |TTT|. We discuss here the stability of the ratios £ = ^pep/t^pp + $pep) 
and rj = $n/($n + $o) among the non-standard solar models which we have built in section 



IV. 

Figs. |29] and |30] show the situations as some input parameter (S pp , Z/X, . . . ) is varied by 
a scaling factor X/X*. In all the presented models, the central temperatures differ by that 
of the starting solar model by no more than 5%. The following features are to be remarked: 

i) variations of £ with respect to the starting model prediction do not exceed ten per 
cent; 

ii) rj is stable to within 20% level or better; 

hi) as in all models we considered the central temperature is smaller or equal to that of 
the starting solar model, the rate of the key reaction p+ 14 N— * 7+ 15 is reduced and the CN 
chain is less equilibrated, resulting in a smaller 15 O neutrino production, which accounts for 
rj > rj* systematically. On the other hand the behaviour of £ looks erratical: it may increase 
as well as decrease when the central temperature is reduced, presumably since £ is sensitive 
to the electron density more than to the temperature. 
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TABLES 



TABLE I. Comparison among several solar models, all calculated including diffusion of helium 
and heavy elements. The "best model with diffusion" of 15] will be used as the Reference Solar 
Model (RSM) in this paper. FRANEC96 indicates our best model with diffusion, model (e) of 
Appendix A. The dag (t) indicates values of c& calculated by us, assuming fully ionized perfect gas 



EOS. The last row indicates the consistency with helioseismology, see section IB 





CGK89 


P94 


DS96 


RVCD96 


RSM 


1 RANEC96 


Ref. 


@ 





@ 


n 


rc Hi 






t©[Gyr] 


4.54 


4.60 


4.57 


4.60 


4.57 


4.57 


L |1U erg/ cm /sj 


O.OZO 


Q QAd 


Q Q/l /I 


O.OOl 


Q Q/l A 


Q Q/l A 


Rq [1U cmj 


n nrnn 

6.9599 


0.9599 


6.960 


6.959 


n nrnn 

6.9599 


6.960 


\ L I -^Iphoto 


U.UZ404 


u.uzoy4 


U.UZZDo 


U.UZUO 


U.UZ440 


n n9/i ^ 


"V" 

-"-in 


0.691 


0.6984 


0. f 295 


0. f012 


0. (024/ 


0.711 


Y 


0.289 


0.2803 


0.2509 


0.2793 


0.27753 


0.269 




0.02 


0.02127 


0.01833 


0.0195 


0.02 


0.0198 


^photo 


0.7265 


0.7290 


0.7512 


0.7226 


0.73507 


0.744 


Y photo 


0.2556 


0.2514 


0.2308 


0.2584 


0.24695 


0.238 


^photo 


0.0179 


0.01964 


0.0170 


0.0190 


0.01798 


0.0182 


Rb/R Q 


0.721 


0.7115 


0.7301 


0.716 


0.712 


0.716 


T b [W 6 K] 


2.142 




2.105 


2.175 


2.204 


2.17 


q,[10 7 cm s _1 ] 


2.21+ 




2.21+ 


2.22+ 


2.25+ 


2.22 


T c [10 7 K] 


1.573 


1.581 


1.561 


1.567 


1.5843 


1.569 


p c [100 gr cm~3] 


1.633 


1.559 


1.554 


1.545 


1.562 


1.518 


Helioseismology 


No 


Yes 


No 


Yes 


Yes 


Yes 



The correspondence 
CEFF = ESai with 



TABLE II. Physical and chemical inputs of the solar models in Table 
between acronyms and references is as follows: C&S70= [156|, MHD = |15 
the Coulomb correction added (see gap ), EFF = Q, BP92 = @, G&N93 = @, G91 = (7|, 
F75 



p4fl, C&F88 = pi. Note that RVCD96 also includes rotational mixing. 



CGK89 P94 



Ref. 



DS96 

m 



RVCD96 
22| 



RSM 



FRANEC96 



OPACITY C&S70 

EOS MHD 

MIXTURE G&N93 
CROSS SECTIONS F75 



OPAL OPAL 

CEFF DS96 

G91 G&N93 

BP92 DS96 



OPAL 

MHD 
G&N93 
C&F88 Table KXVIII 



OPAL 
BP92 
G&N93 



OPAL 
OPAL 

G&N93 



Table KXVIII 
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TABLE III. Comparison among the neutrino fluxes of the solar models of Table |. All of them, 
except for DS96, are SSMs according to our definition. 



Ref. 



P94 



DS96 



RVCD96 



RSM 



FRANEC96 



$ pp [10 9 cm" 2 s" 1 ] 59.1 61.0 59.4 59.1 59.92 

$pe P [10 9 cm" 2 s- 1 ] 0.139 0.143 0.138 0.140 0.14 

$ Be [10 9 cm" 2 s- 1 ] 5.18 3.71 4.8 5.15 4.49 

$ N [10 9 cm" 2 s- 1 ] 0.64 0.382 0.559 0.618 0.53 

$ [10 9 cm~ 2 s" 1 ] 0.557 0.374 0.481 0.545 0.45 

$ B [10 6 cm" 2 s" 1 ] 6.48 2.49 6.33 6.62 5.16 



S Ga [SNU] 
Scz[SNU] 



136.9 
9.02 



115 
4.1 



132.77 
8.49 



137.0 
9.3 



128 
7.4 



TABLE IV. Helioseismological determinations of the present surface He abundance, Y p h t - In 



RVCD96 uncertainties depending 


on the inversion method and on the EOS are included. 


Reference 




^ photo 


Dappen (1988) p0| 




Y=0.233 ± 0.003 


Dappen (1991) pH 




Y=0.268 ± 0.002 


Dziembowski (1991) | 


155] 


Y=0.234 ± 0.005 


Vorontsov (1991) |62 


t 


Y=0.250 ± 0.010 


Dziembowski (1994) | 


w 


Y=0.24295 ± 0.0005 


Hernandez (1994) |6 


1 


Y=0.242 ± 0.003 


RVCD96 HI 




Y=0.250 ± 0.005 



TABLE V. For each component of the neutrino flux, we show the average neutrino energy {E v ) 
and the averaged neutrino capture cross sections (10 -9 SNU cm 2 s = 10 -45 cm 2 ) for Chlorine and 
Gallium, with errors at la level . When averaging the pp and pep components to get p, we use 
the relative weights of the RSM; similarly for 13 N and 15 to get CNO. All data from [24| but for 



a C l,B from (165 1 









&Ga,i 




[MeV] 


[10" 9 SNU cm 2 s] 


[10" 9 SNU cm 2 s] 


PP 


0.265 


0. 


1.18(1 ±0.02) 


pep 


1.442 


1.6 (1±0.02) 


21.5 (1±0.07) 


p=pp+pep 


0.268 


0.0038(1 ± 0.02) 


1.23(1 ±0.02) 


7 Be 


0.814 


0.24 (1±0.02) 


7.32(1 ±0.03) 


13 N 


0.707 


0.17 (1±0.02) 


6.18(1 ±0.03) 


15q 


0.996 


0.68 (1 ± 0.02) 


11.6 (1±0.06) 


CNO= 13 N + 15 


0.842 


0.41 (1±0.02) 


8.72(1 ±0.05) 


8 B 


6.71 


1110. (1±0.03) 


2430. (1 ±0.25) 
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TABLE VI. The main characteristic of each neutrino experiment: type, detection reaction, 
energy threshold E t h, experimental results with statistical and systematical errors. In the last 



column the Reference Solar Model [15| predictions are presented. Errors are at la level. 



Experiment 


type 


TP a 


result b 


BP95 b 


Homestake 


radiochemical 
u+ 37 Cl^ e~+ 37 Ar 


0.814 


2.55 ±0.17 ±0.18 




KAMIOKANDE 


OvfCbli UC1 .J. 1 1 .£-, 

v + e~ — > v + e~ 


7 


(2 73 ± 17 ± 341 


6 62(1 QQ + %-}$) 


GALLEX 


radiochemical 
i/+ 71 Ga^ e"+ 71 Ge 


0.233 


77.1 ± j 


137t* 


SAGE 


radiochemical 
j/+ 71 Ga^ e"+ 71 Ge 


0.233 


69 ± 10±f 


137t 8 7 



a Energy in Mev 

b in SNU for radiochemical experiments; in 10 6 cm _2 s _1 for KAMIOKANDE. 



TABLE VII. Contribution from the main components of neutrino flux to the signals (SNU) in 
71 Ga and 37 C1 detectors according to the RSM, from [15|. 



71 



Ga 



37 



CI 



pp 


69.7 


0.0 


pep 


3.0 


0.22 


7 Be 


37.7 


1.24 


13 N 


3.8 


0.11 


15Q 


6.3 


0.37 


8 B 


16.1 


7.36 



Total 



136.6 



9.30 



TABLE VIII. Experimental information on the fluxes of intermediate energy neutrinos (units 
of 10 9 cm _2 s _1 ). In (a) only the luminosity constraint is assumed. In (b) we assume also £ and rj 



as given by the RSM |15[ , see Sec. IIP . The best fit points and upper limits, within 2a from each 
experimental results, are presented. 







(a) 






(b) 






13 N 


7 Be 


15Q 


pep 


7 Be±CNO 


CNO 


best fit 


-2.2 


-1.6 


-0.8 


-0.2 


-2.2 


-1.7 


upper limit (2a) 


1.8 


1.5 


0.8 


0.5 


0.9 


0.7 
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TABLE IX. Information on neutrino fluxes. $g e and ^CNO (^b) ^ n un its °f 10 9 cm 2 s 



(10 6 cnT 



All bounds are at the 99.5% C.L. Direct information (exp) is only available for 



B neutrinos, from Ref. M, the indicated error here corresponds to 3a. The bounds in (bl) 
correspond to no prior knowledge on the unknown variables. In (b2) we assume a priori $j > 0. 
The results of SSMs are also shown (same notation as in Table |). 



Ref. 



(exp) 
§ 



(bl) 



(b2) 



P94 



RVCD96 



RSM 



FRANEC96 



$B 

$ Bc+CNO 



2.73±1.14 



Be 
CNO 



<0.7 
<0.7 
<0.5 



<2.0 
<2.0 
<1.5 



6.48 
6.38 
5.18 
1.20 



6.33 
5.9 
4.8 
6.33 



6.62 
6.31 
5.15 
1.16 



5.16 
5.47 
4.49 
0.98 



TABLE X. Results of variations of S pp , for quantities characterizing the solar interior and for 
the neutrino fluxes. We present the calculated values of a, see Eqs. flS9| ) and (p6[), for the case 
of small and large variations. For this latter case also the variances Aa are shown, Eq. (^). For 
the dependence on the central temperature, the corresponding (3 coefficients, Eq. (|103j ), are also 
presented. 





0.9-1.1 


1-3.5 


0.9-1.1 


1-3.5 


coefficient 


a 


a ± Aa 





P± A/? 


T 


-0.11 


-0.11 ± 0.001 






P 


-0.37 


-0.37 ± 0.01 


3.3 


3.3 ± 0.05 


P 


-0.49 


-0.51 ± 0.01 


4.4 


4.6 ± 0.1 


R 


+0.12 


+0.13 ± 0.003 


-1.1 


-1.1 ± 0.03 


X 


-0.01 


-0.03 ± 0.004 


0.1 


0.3 ± 0.04 


$p 


+0.11 


+0.07 ± 0.01 


-0.9 


-0.6 ± 0.1 


$Be 


-1.02 


- 1.1 ± 0.04 


8 


9 ± 0.3 


$ CNO 


-2.7 


-2.2 ± 0.3 


21 


18 + 2 


$B 


-2.7 


-2.7 ± 0.1 


21 


22 ± 1 
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TABLE XL Variations of opacity. Same notation as in Table |X[ 



opa 


0.9-1.1 


0.6-1 


0.9-1.1 


0.6-1 


coefficient 


Q: 


a ± Aa 




P±A0 


T 


+0.12 


+0.12 ± 0.003 


- 


- 


P 


-0.06 


-0.05 ± 0.006 


-0.5 


-0.4 ± 0.05 


P 


-0.08 


-0.07 ± 0.003 


-0.7 


-0.6 ± 0.03 


R 





+0.01 ± 0.01 





+0.1 ± 0.1 


X 


-0.28 


-0.28 ± 0.01 


-2.3 


-2.3 ± 0.1 


$p 


-0.11 


-0.09 ± 0.01 


-0.9 


-0.68 ± 0.1 


$Be 


+1.1 


+1.2 ± 0.1 


+8.5 


+9.1 ± 1 


^CNO 


+1.7 


+1.7 ± 0.1 


1 3 4 


-1-13 + 1 




+2.4 


+2.6 ± 0.2 


19 


20.1 ± 1 




TABLE XII. Variations of Z/X. Same notation 


as m lable [A] 




(Z/X)/(Z/X)* 


0.9-1.1 


0.1-1 


0.9-1.1 


0.1-1 


coefficient 


a 


a ± Aa 




P± A/3 


T 


+0.06 


0.05 ± 0.005 


- 


- 


P 


+0.03 


+0.02 ± 0.003 


0.5 


0.4 ± 0.06 


P 


0.006 


+0.001+ 0.001 


0.1 


0.02 ± 0.02 


R 


-0.04 


-0.01 ± 0.004 


-0.7 


-0.2 ± 0.1 


X 


-0.2 


-0.13 ± 0.01 


-3 


-2.6 ± 0.2 


$p 


-0.06 


-0.04 ± 0.01 


-0.9 


-0.7 ± 0.1 


$ Be 


+0.62 


+0.54 ± 0.03 


9.9 


10.6 ± 0.4 


$ CNO 


+2 


+1.7 ± 0.1 


31 


33 + 2 




+1.3 


+1.1 ± 0.1 


21 


21 ± 1 


TABLE XIII. Variations of solar age. Same notation as in Table |X[ 


t©/t0 


0.9-1.1 


0.1-1 


0.9-1.1 


0.1-1 


coefficient 


a 


a ± Aa 


(3 


P±Ap 


T 


+0.03 


+0.02 ± 0.007 






P 


+0.18 


+0.11 ± 0.04 


+6 


+5.5 ± 2 


P 


+0.12 


+0.08 ± 0.03 


+4 


+4 + 2 


R 


-0.08 


-0.05 ± 0.01 


-3 


-2.5 ± 0.5 


X 


-0.2 


-0.13 ± 0.05 


-7 


-6.5 ± 2.5 


$p 


-0.08 


-0.04 ± 0.01 


-1.4 


-0.8 ± 0.1 


$Be 


+0.57 


+0.5 ± 0.1 


+10 


+11 ± 1 




+0.9 


+0.5 ± 0.2 


+16 


+12 ± 3 


$B 


+1. 


+0.8 ± 0.2 


+18 


+20 ± 1 
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TABLE XIV. The power law coefficients for the reaction rates as a function of temperature, 
Eq. (glf), calculated for T = 15.6 • 10 6 K. 



■ ■ Alnlav) 

reaction 7 = gr^r 

P+P 4 

3 He + 3 He 16 

3 He + 4 He 16 

p+ 7 Be 13 

e+ 7 Be -0.5 

p+ 14 N 20 



TABLE XV. The (3 coefficients connecting the neutrino fluxes with the temperature, Eq. (103). 
The components of neutrino flux are indicated in the first column. The values presented are the 
best fit to the numerical calculations performed when each input parameter is varied in the range 
specified in the second row. 



Parameter 


u pp 


Opacity 


Z/X 


t 


Scaling factor 


1 - 3.5 


0.6 - 1 


0.1 - 1 


0.1 - 1 


$PP 


-0.6 


-0.7 


-0.7 


-0.8 


^pep 


2.2 


-2.3 


-1.7 


0.5 


^P = ^pp+pep 


-0.6 


-0.7 


-0.7 


-0.8 


$Bc 


9 


9 


11 


11 




15 


12 


31 


9 




24 


15 


36 


18 


$CNO = $N + $0 


18 


13 


33 


12 




22 


21 


20 


20 



TABLE XVI. Values of the coefficients olt,x relating temperature to the input parameter, 
T c = T*(X/X*) aT ' x see Eq. (103), our estimated uncertainties of the input parameters (AX/X), 
and variations (5X/X) required to reduce T c by 3%, 7% and 13%, respectively. 



X 




AX 

X 


sx 

X 


-3% 


SX 
X 


-7% 


sx 
x 


-13% 


Spp 


-0.13 


1% 


+25% 


+70% 


+190% 


Opacity 


+0.13 


5% 


-20% 


-40% 


-66% 


Z/X 


+0.06 


10% 


-40% 


-70% 


-90% 


T 


+0.05 


0.6% 


-45% 


-75% 


-94% 
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TABLE XVII. Estimates for S33 (0) [MeVb]. 

CF88 5.57 

PA91 @ 5.0+0.3 

All data 5.3+0.2 

£ > lOOKeV 5.2+0.2 

adiabatic scr. 5.1+0.2 



TABLE XVIII. Estimates for S* 34 (0) [keVb]. 



CF88 || 
PA91 H 
quadratic 
exponential 



0.54 
0.533+0.017 
0.48+0.01 
0.51+0.01 



TABLE XIX. Variation of x = 5 , 34/V5'33- We present the coefficients ctj for the parameteriza- 
tion <3?j = f° r X i n the range indicated in the first row. In the second column the values 
of Ref. |24j] are shown. For the case of large variations, the variances ±Aa are also presented. 



x/x* 


0.9-1.1 


0.9-1.1 


0.1-1 


reference 


@ 


this work 


this work 


$p 


-0.06 


-0.05 


-0.05 ± 0.01 


$Be 


+0.86 


+0.86 


+0.92 ± 0.02 


$ CNO 


-0.05 


-0.04 


-0.02 ± 0.01 


$B 


+0.80 


+0.92 


0.91 ± 0.02 



TABLE XX. Different determinations of 5(0)i 7 , from jl20t| 



Ref. 
[p5 | 

mi 
m\ 

M\ 

El9| 



5(0)i7 [eV b] 
15 + 6 
27 + 4 
25.2 + 2.4 

19.4 + 2.8 

41.5 + 9.3 
20.2 + 2.4 
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TABLE XXI. Power law coefficients for variations of S17. Same notation as in Table 



S 17 /S{ 7 0.1-10 

¥^ (+1.1 ±0.4)10" 4 

$ Bc (-2.5±2.0)10~ 3 
<J? CN0 (-1.7±0.7)10- 4 
$ B (0.996 ± 0.003) 



TABLE XXII. Comparison among solar models with different screening predictions: NOS=no 



screening, WES=weak screening |3§|, MIT=Mitler 1977 [140] , GDGC=Graboske et al. 1973 |3g | 
and CSK=Carraro e£ aZ. 1988 [139]. We show the central temperature T C [10 7 K], the helium 



abundance in mass Y, the metal fraction Z, the values of each component of the neutrino flux 
[10 9 cm~ 2 s _1 ], the calculated signals for the Chlorine (CI) and the Gallium (Ga) experiments [SNU], 
from [P0l. 





NOS 


WES 


MIT 


GDGC 


CSK 


T c 


1.573 


1.566 


1.566 


1.564 


1.567 


Y 


0.288 


0.289 


0.289 


0.289 


0.289 


Z (xlO 2 ) 


1.85 


1.85 


1.84 


1.84 


1.85 


pp 


60.0 


59.6 


59.7 


60.0 


59.7 


pep 


0.146 


0.142 


0.142 


0.143 


0.143 


7 Be 


4.82 


4.97 


4.93 


4.79 


4.94 


8 B (xlO 3 ) 


5.51 


6.36 


6.13 


5.59 


6.21 


13 N 


0.46 


0.55 


0.52 


0.47 


0.54 


15q 


0.39 


0.48 


0.45 


0.40 


0.47 


CI 


7.7 


8.8 


8.5 


7.8 


8.6 


Ga 


130 


134 


133 


130 


134 



TABLE XXIII. Enhancement factors for different screening prescriptions, calculated in central 



solar conditions. Same notations of Table XXII, from 



WES 



MIT 



GDGC 



CSK 



p + p 
He+He 
7 Be+p 
14 N+p 



1.049 
1.213 
1.213 
1.403 



1.045 
1.176 
1.171 
1.293 



1.049 
1.115 
1.112 
1.192 



1.038 
1.158 
1.169 
1.324 



85 



TABLE XXIV. Best fits to the combined experimental results obtained by varying the seven 
possible combinations of the parameters T c , S^ s and S17. For comparison, the first two rows also 
show the experimental and RSM results. The first column reports the x 2 /d.o.f (just the \ 2 when 
the number of parameters is greater or equal to the number of data); a dag (f) indicates that the 
value is not a local minimum, but it is the lowest value within the explored region (Sn < 5Syf M 
and 633 s < 200S§f M ). The second column reports the best-fit values of the parameters in units 
of the RSM, and the last three columns the corresponding signals for the Gallium, Chlorine and 
KAMIOKANDE experiments. 





X 2 /d.o.f. 


best-fit 
values 


Saa 


[SNU] 


Sci 


10 6 ernes' 1 


experiment 






74 ±8 


2. 


55 ±0.25 


2.73 ± 0.38 


RSM 


58.5/3 




137.0 




9.3 


6.62 


Sn 


50.7/2 


0.28 


125 




4.00 


1.85 


T c 


22.0/2 


0.936 


103 




2.71 


1.54 


ores 
°33 


18.1/2 


13.9 


99 




2.90 


1.72 


T c 

ores 
^33 


16.1/1 


0.973 
4 


97 




2.68 


1.66 


ores 
^33 

Sl7 


13.9t/l 


200 
4.2 


96 




3.03 


1.96 


T c 
Sn 


10.7T/1 


0.879 
5.0 


94 




2.74 


1.94 


T c 

ores 
^33 

Sn 


7.4t 


0.929 
9 

5.0 


90 




2.76 


2.07 
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TABLE XXV. The best fits when just two experimental results are included. Same notations 
as in Table XXIV. 





Ga 


+ Ka 




CI 4- 


Ga 




CI + Ka 




X /d.o.t. 




best 


X /d.o.t 


best 


X 2 /d.o.f. best 


V> Ql\/T 
rtoM 


on c /O 

oy.o/ Z 






oU.z/ z 






41.9/2 


o 

<->17 


35.0/1 




0.50 


30.5/1 


0.12 




22 8/1 21 


T c 
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TABLE XXVI. Comparison among 


solar models 


obtained with different versions of the 


FRANEC co 


de. The labels I 


a) to (e) corresponds to the models defined 


in the Appendix A. 


Our best Standard Solar Model is (e). The last column shows the helioseismological results. 






(a) 


(b) 


(c) 


(d) 


(e) 


Helioseism. 














our best 




t©[Gyr] 




4.57 


4.57 


4.57 


4.57 


4.57 




L [1O 33 erg/cm 2 /s] 


3.846 


3.843 


3.844 


3.843 


3.844 




R [1O 10 cm] 




6.961 


6.963 


6.959 


6.959 


6.960 




(Z/X)photo 




0.0245 


0.0245 


0.0245 


0.0245 


0.0245 




a 




2.023 


1.774 


1.786 


1.904 


1.901 




Xj n 




0.699 


0.718 


0.722 


0.711 


0.711 




Y 

L in 




0.284 


0.265 


0.261 


0.269 


0.269 




^in 




0.0171 


0.0176 


0.0177 


0.0198 


0.0198 




^photo 




0.699 


0.718 


0.722 


0.743 


0.744 




i photo 




0.284 


0.265 


0.261 


0.238 


0.238 


0.233 - 0.268 


^photo 




0.0171 


0.0176 


0.0177 


0.0182 


0.0182 




Rb/R Q 




0.738 


0.726 


0.728 


0.716 


0.716 


0.710 - 0.716 


T 6 [10 6 K] 




1.99 


2.10 


2.08 


2.17 


2.17 




Q,[10 7 cm s x ] 


2.11 


2.16 


2.16 


2.22 


2.22 


2.21 - 2.25 


T c [10 7 K] 




1.555 


1.545 


1.542 


1.569 


1.569 




p c [100 gr cm" 


-3] 


1.524 


1.472 


1.470 


1.514 


1.518 




Y c 




0.63 


0.61 


0.61 


0.63 


0.63 
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TABLE XXVII. Neutrino fluxes and signals obtained with different versions of the FRANEC 
code. Labels (a) to (e) correspond to the models defined in Appendix A. Our best prediction is 
(e). 





(a) 


(b) 


(c) 


(d) 


(e) 












our best 


$ pp [10 9 cm" 2 s- 1 ] 


60.17 


60.37 


60.66 


59.76 


59.92 


$ pep [10 9 cm" 2 s" 1 ] 


0.14 


0.14 


0.14 


0.14 


0.14 


$ Bc [10 9 cm- 2 s" 1 ] 


4.58 


4.22 


4.09 


4.71 


4.49 


% [10 9 cm" 2 s- 1 ] 


0.39 


0.36 


0.35 


0.52 


0.53 


$ [10 9 cm" 2 s- 1 ] 


0.33 


0.30 


0.29 


0.45 


0.45 


$ B [10 6 cm~ 2 s- 1 ] 


4.73 


4.18 


3.95 


5.37 


5.16 


S Gcl [SNU] 


126 


121 


120 


130 


128 


S CZ [SNU] 


6.9 


6.2 


5.9 


7.7 


7.4 



TABLE XXVIII. Astrophysical S-factors [MeV barn] and their derivatives with respect to 



energies S' [barn] for the RSM [15] and for our models. 





RSM 


models (a-d) 


model (e) 








our best 


S(0)u 


3.89xl0" 25 


3.89 xlO" 25 


3.89 xlO" 25 


S"(0)u 


4.52xl0- 24 


4.52 xlO" 24 


4.52 xlO" 24 


5(0)33 


4.99 


5.00 


5.1 


S"(0) 33 


-0.9 


-0.9 


3.0 


5(0)34 


5.24 xlO- 4 


5.33 xlO" 4 


5.1 xlO" 4 


5'(0) 34 


-3.1 xlO" 4 


-3.10 xlO" 4 


-4.23 xlO" 4 


5(0) 17 


2.24 xlO" 5 


2.24 xlO" 5 


2.24 xlO" 5 


5'(0) 17 


-3.00 xlO" 5 


-3.00 xlO" 5 


-3.00 xlO" 5 


S{0) 12C + p 


1.45 xlO- 3 


1.40 xlO" 3 


1.40 xlO" 3 


s'(o) 12C + p 


2.45 xlO" 4 


4.24 xlO" 3 


4.24 xlO" 3 


S{0) 13C + p 


5.50 xlO- 3 


5.50 xlO" 3 


5.50 xlO" 3 


S'(0) 13C + p 


1.34 xlO" 2 


1.34xl0- 2 


1.34xl0- 2 


5(0)i4_/y + p 


3.29 xlO" 3 


3.32xl0- 3 


3.32xl0- 3 


5 / (0)i4jv + p 


-5.91 xlO- 3 


-5.91xl0- 3 


-5.91xl0- 3 


S ( )i6N(p^) 16 O 


6.40xl0- 2 


6.40xl0- 2 


6.40xl0- 2 


S '(°)i S N(p, 7 ) 16 


3.00xl0- 2 


3.00 xlO" 2 


3.00 xlO" 2 


S ^^N(p,a) 12 C 


7.80x10 


7.04 xlO 


7.04 xlO 


S '^N(p,a) 12 C 


3.51xl0 2 


4.21xl0 2 


4.21xl0 2 


5(0)i 6 Q _|_ p 


9.40xl0- 3 


9.40xl0- 3 


9.40xl0- 3 


S'(0) W Q + p 


-2.30xlQ- 2 


-2.30xlQ- 2 


-2.30xl0- 2 
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FIGURES 



FIG. 1. The pp chain. The probability of the different branches are from the Reference Solar 
Model [15]. The neutrino energies E„ are also indicated. 



FIG. 2. The CNO cycle. 



FIG. 3. The photospheric helium mass fraction photo an d the depth of the convective zone 
(R*/Ro): 



a) as constrained by helioseismology (the dotted rectangle), see Sect. [B; 

b) as predicted by solar models without diffusion, open circles from top to bottom correspond to 
@,^T],||Jl^,|T|,|l5|,|§; 

c) as predicted by solar model with helium diffusion, full squares from top to bottom correspond 
to 000); 

d) as predicted by solar model with helium and heavy elements diffusion: the full circles, from top 



to bottom correspond to [22,45, 16, 56, 23 1) and the full diamond indicates the RSM pi 



FIG. 4. The solar neutrino spectrum, from []2J]. For continuous sources, the differential flux is 
in cm~ 2 s _1 MeV -1 . For the lines, the total flux is in cm -2 s _1 . 



FIG. 5. For the indicated components, df is the fraction of neutrinos produced inside the sun 
within dR. On the bottom (top) scale the radial (mass) coordinate is indicated. 

FIG. 6. The 8 B and 7 Be neutrino fluxes, consistent with the luminosity constraint and exper- 
imental results, for standard neutrinos. The dashed (solid) lines correspond to central (±lcx) ex- 
perimental values for Chlorine (CI), Gallium (Ga) and KAMIOKANDE (Ka), see Eqs. @|26||2|). 
The hatched area corresponds to the region within 2a from each experimental result. The di- 
amond represents the prediction of the Reference Solar Model |f[q| , and the bars the estimated 
uncertainties. 



FIG. 7. Same as in Fig. ^, for the 8 B and 7 Be + CNO neutrino fluxes. The luminosity 
constraint is supplemented with the estimates for < 5pep/( < ^pep + ^pp) an d <J?n/(<J>n + ^o)' f rom 



the RSM [15]. The prediction of the RSM model |L5| (full diamond) is shown, together with those 
of other solar models including diffusion of helium and heavy elements (full circles), corresponding 

Solar model calculations without diffusion are represented by 



from right to left to |T§g|6],[S 
open circles, corresponding to [JO* IS ,2C,17|, again from right to left. 



FIG. 8. The fate of 7 Be nuclei. 
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FIG. 9. For x = S pp /S* p = 2 the behaviour of several structure parameters, Oj(x, m)/f2j(m)*, 
as a function of the mass coordinate in the whole internal radiative region. The considered struc- 
tural parameters are: radius (R), density (p), temperature (T), pressure (P) and the hydrogen 
mass fraction (X). 

FIG. 10. Same as in Fig. ||, for opa = 0.7. 

FIG. 11. Same as in Fig. g for (Z/X)/(Z/X) 5SM =0.5. 

FIG. 12. The temperature profiles T(m) normalized to T(m)* for models with the indicated 
values of the solar age. 

FIG. 13. Same as in Fig. ^, for the model with t© = 4 Gyr. 

FIG. 14. The temperature profiles T(m) normalized to T SSM (m) for a few representative 
non-standard solar models, from |7J 



FIG. 15. The behaviour of ^pp, anc ^ as a function of the central temperature T c when 
varying S pp , opacity, Z/X and age, from ]73[ . 



FIG. 16. Relations among the temperatures T at the 7 Be and pep peak production zones 
(R/R Q = 0.06 and R/R = 0.09, respectively) and the central temperature T c in non-standard 



solar models. Data from numerical calculations are shown with the same symbols as in Fig. 15, 



while full lines show the homology relations T{ = T c (T*/T*), from [73]. 



FIG. 17. The 3 He ( 3 He , 4 He )2p reaction. Experimental data for the S-factor as a function of 
CM energy (lower scale). The corresponding temperature Tg(= T ■ 10~ 9 ), such that E (Tg) = E cm , 



is also indicated (upper scale). The full curve corresponds to Eq. (132). The Gamow peak for the 
solar center is indicated by the arrows. The listed symbols corresponds to data in |9l|-p6|, from 
top to bottom. 

FIG. 18. The 3 He ( 4 He ,7) Be reaction. Experimental data for the S-factor as a function of 
CM energy (lower scale). The corresponding temperature T9, such that E (Tg) = E cm , is also 
indicated (upper scale). The full (dashed) curve corresponds to Eq. (133) (134)- The listed 



symbols corresponds to data in | 105| , p9f - [l0llJl03 , 106 , |102 , 10^ , 107], from top to bottom. 



FIG. 19. The central temperature T c as a function of 533. 
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FIG. 20. The flux of (a) 7 Be neutrinos and (b) 8 B neutrinos as a function of % = Szi/ Ss-s' f° r 
several values of S33, as calculated by using the FRANEC code. The dashed line corresponds to 
Eq. flgg), from pi. 



FIG. 21. Proton capture (a) and Coulomb dissociation (b) reactions. 



FIG. 22. The central temperature T c as a function of S17. 



FIG. 23. Enhancement factors along the solar profile. The results of weak screening |136| ] 
(dashed curves), Graboske et al. 1973 [ |138| (dot-dashed curves), Mitler 1977 |140[| (solid curves) 
and Carraro et al. 1988 |13S| 1 (dotted curves) are shown, from 



FIG. 24. Dependence of the main solar fluxes on T c , Sn and 6*33 (this latter parameterized 
by its zero-energy resonant contribution Sjjf 5 ). The coefficients Qi are defined in Eq. 



FIG. 25. 7 Be plus CNO neutrino fluxes vs. 8 B neutrino flux. The three stripes, labeled Ga, 
CI and Ka, confine the regions allowed at the 2a level by the three current experimental data 
(GALLEX + SAGE, Chlorine, and KAMIOKANDE) with the constraint due to the luminosity 
sum rule. The hatched area emphasizes the interception of the three regions. The diamond shows 
the RSM prediction. The solid (dashed, dotted) line shows the effect of decreasing T c (increasing the 
resonant part of S33, decreasing S17). Dots indicate the fluxes at specific values of the parameters 
shown by the label in units of their RSM values. 



FIG. 26. Similar to Fig. 25, but this time two parameters are changed in the same model. The 
temperature is decreased down to T c = 0.88 (solid curve) and, then, Sn is increased up to 5 times 
its RSM value (dotted line). 



FIG. 27. Similar to Fig. 25, but this time all three parameters are changed in the same model. 
First S^ s is increased from zero up to 9 times the RSM value of S33 (dashed curve), then T c is 
decreased down to 0.929 (solid curve) and, finally, Sn is increased up to 5 times its RSM value 
(dotted line). The order of the transformations is obviously inconsequential. 



FIG. 28. The dotted area shows the result of arbitrarily varying the enhancement factors ip+p, 
/He+He an d /p_|_i4N up to / = 6. We recall that in standard solar model calculations / ^ 1.2, see 
sect. 0. 

FIG. 29. The behaviour of £ = < &pep/( < ^pp + *£pep) as one of the input parameters is varied by 
a scaling factor X/X*. 
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FIG. 30. The behaviour of 77 = 3>n/(<I>n + 3>o) as one of the input parameters is varied by a 
scaling factor X/X*. 
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